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ABSTRACT 


This paper presents a review of studies on mathematical modeling of solid oxide fuel cells (SOFCs) with 
respect to the tubular and planar configurations. In this work, both configurations are divided into five 
subsystems and the factors such as mass/energy/momentum transfer, diffusion through porous media, 
electrochemical reactions with and without CO oxidation, shift and reforming reactions, and polarization 
losses inside the subsystems are discussed. Using variety of fuels fed to SOFCs is issued and their effect 
on the system is compared briefly. A short review of solid oxide fuel cell configurations and different 
flow manifolding are also presented in this study. Novel models based on statistical data-driven approach 
existing in the literatures are considered shortly. Although many studies on solid oxide fuel cells modeling 
have been done, still more research needs to be done to improve the models in order to predict the fuel cell 
behaviors more accurately. At the end of this paper the works and studies that can be done for improving 


Polarization the fuel cell models is suggested and pointed by the authors. 


Novel SOFC models © 2011 Elsevier Ltd. All rights reserved. 
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1. Introduction 


Ever increasing energy consumption, rising public awareness 
for environmental protection and higher prices of fossil fuels have 
motivated many to look for alternative/renewable energy sources. 
Fossil fuel resources are limited and are expected to end early in 
the next century. By 2015, the world fossil fluid fuel demand will 
exceed the world fossil fluid fuel production, which is expected to 
lead to an energy shortage crisis unless a sustainable alternative 
fuel will be available by them [1-3]. 

Small-scale power-generation systems such as wind turbines, 
photovoltaics, micro-turbines and fuel cells can play an important 
role in meeting consumer demands for “greener” energy. These 
small-scale generation systems, usually located near customers 
rather than at central or remote locations, can be used to form 
distributed power generation networks. Among different types of 
small-scale power generation systems, fuel cells have received 
more attention, because they can provide both heat and power. 
Household fuel cell lets households to generate part of the electric- 
ity they require, reducing the electricity purchase from an electric 
power company. In addition, they serve as a cogeneration system, 
which use thermal energy produced by electricity generation for 
water and air heating. It would be best if exhaust heat emitted dur- 
ing the generation process from large-scale power stations could 
also be converted to thermal energy for heating air or water to use 
the inputted energy more efficiently. However, it is not considered 
practical to send hot water (thermal energy) to homes by pipelines 
from far-removed locations of such large-scale power stations. 

Fuel cells are static energy conversion devices that partially con- 
vert the chemical energy of fuels directly into electrical energy and 
produces water as its byproduct [4,5]. The basic principle of the fuel 
cell was discovered in the year 1838 by Swiss scientist Christian 
Friedrich Schnbein. In 1839 Sir William Grove developed the first 
fuel cell based on reversing the electrolysis of water by an accident 
[3]. In 1950 Francis Bacon at Cambridge University demonstrated 
the first 5kW alkalinahoe fuel cell. After the successful develop- 
ment of alkaline fuel cells, NASA needed a compact system to 
generate electricity for space shuttle applications. In 1970s, inter- 
national fuel cells developed a 12 kW alkaline fuel cell for NASA’s 
space shuttle orbiter to supply reliable power without the use of 
any backup powers like batteries. Beginning in the mid-1960s, the 
research work was focused on further development of various fuel 
cells for applications like stationary powers and transportations. 
Further the government agencies in the USA, Canada, Japan and 
Malaysia have significantly increased their funding for fuel cell in 
R&D [6]. 

There are six main types of fuel cell: alkaline fuel cell (AFC), 
phosphoric acid fuel cell (PAFC), molten carbonate fuel cell 
(MCFC), polymer electrolyte membrane fuel cell (PEMFC) and direct 
methanol fuel cell (DMFC), solid oxide fuel cell (SOFC). Table 1 
depicts their characteristics briefly. 

After the development of power conversion devices, much more 
researches are doing on solid oxide fuel cells to reduce its higher 
installation cost. Solid oxide fuel cells are considered to be one of the 
most advanced designs for mid- to large-scale applications. They 
are among the most promising types of fuel cells currently being 
considered as a power source for automobiles and stationary power 


plants [2,3,7]. Because the electrolyte is a layer of ceramic material 
with high-temperature endurable porous-media electrodes, SOFCs 
can generally operate in a high-temperature range (800-1000 C). 
High operating temperature has some advantages, such as high 
energy conversion efficiency, flexibility of usable fuel type, and high 
temperature exhaust gas. Disadvantages include potential ther- 
mal fatigue failure of the cell material and sealing under the high 
temperature, as well as the fact that cell temperature fluctuations 
induce thermal stress in the cell ceramics. Thus, it is important to 
operate SOFCs in such a way that the stack temperature remains 
within a tight design range. 

In this case, modeling is a low cost method for studying and 
investigating the fuel cell in order to optimize and control the fuel 
cell behavior, increase its efficiency and performance, and finally 
decrease the higher installation cost. Among many mathematical 
modeling studies, there are some review papers, concerning differ- 
ent aspects of solid oxide fuel cells investigation which introduce a 
complete view of all the proceeding works to the researchers who 
are going to develop and improve them. In this regard, Badwal et al. 
[8] reviewed the different kinds of materials used in electrolyte and 
electrodes. Similar work has been done by Zhu et al. [9], presented a 
review paper about anode materials used in SOFCs. The criteria for 
the anode of SOFC were first discussed. They also discussed about 
the problems existing in the developed anode materials. In this 
paper the advantages and disadvantages of anodes materials were 
illustrated in order to give a view for further research and develop- 
ment of new generation of anode materials for SOFC. A. Stambouli 
[6] reconsidered SOFCs as environmentally clean and an efficient 
source of energy. The advantages and disadvantages of models 
regarding to AC impedance was studied completely with Q. Huang 
et al. [10]. Zhang et al. [11] performed a literature study about dif- 
ferent strategies and concepts for SOFC-based integration systems, 
which are timely transformational energy-related technologies 
available to overcome the threats posed by climate change and 
energy security. A complete review paper about dynamic models of 
single-cell SOFCs with respect to various methods, assumptions and 
models in design, diagnosis and SOFC operation was presented by 
Bhattacharyya et al. [12] but it was less discussion on steady-state 
models. The impact of the cell geometry was also presented in their 
review paper. S. Kakaç et al. [13] reconsidered the numerical models 
of SOFCs and summarized the present status of the SOFC modeling. 
V. Lawlor et al. [14] discussed in two parts of a review paper about 
micro-tubular (MT-SOFC) solid oxide fuel cell. Part I illustrated the 
reader to the MT-SOFC stack and its applications, indicating who is 
researching what in this field and also specifically investigate the 
design issues related to multi-cell reactor systems called stacks. 
Part II [15] reviewed in detail the combinations of materials and 
methods used to produce the electrodes and electrolytes of MT- 
SOFC’s. Also the role of modeling and validation techniques used 
in the design and improvement of the electrodes and electrolytes 
was investigated. M. Andersson et al. [16] presented a literature 
study with respect to various transport processes such as mass, 
heat, momentum, and also charge, electrochemical, shift and inter- 
nal reforming reactions. This work focused on SOFCs designing by 
combining the accuracy at microscale with the calculation speed at 
macroscale, based on clear understanding of transport phenomena, 
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Table 1 
Types of fuel cell [1]. 


Fuel cell Temperature (°C) Efficiency (%) Application Advantage Disadvantages 
Alkaline fuel cell (AFC) 50-90 50-70 Space application High efficiency Intolerant to CO2 in impure H2 
and air, corrosion, expensive 
Phosphoric acid fuel cell (PAFC) 175-220 40-45 Stand-alone & combined Tolerant to impure H2 Low power density, corrosion 
heat & power & sulfur poisoning 
Molten carbonate fuel cell (MCFC) 600-650 50-60 Central, stand-alone & High efficiency Electrolyte instability, 
combined heat & power corrosion & sulfur poisoning 
Solid oxide fuel cell (SOFC) 800-1000 50-60 Central, stand-alone & High efficiency & direct High temperature, thermal 
combined heat & power fossil fuel stress failure, coking & sulfur 
poisoning 
Polymer electrolyte membrane 60-100 40-50 Vehicle High power density, low Intolerant to CO in impure H2 
fuel cell (PEMFC) temperature and expensive 
Direct methanol fuel cell (DMFC) 50-120 25-40 No reforming, high power Low efficiency, methanol 


density & low temperature crossover & poisonous 


byproduct 


chemical reactions and functional requirements. Useful methods 
were investigated for SOFCs modeling. Coupling methods between 
different approaches and length scales by multiscale models were 
outlined in order to enhance the understanding for detailed trans- 
port phenomena and making a correct decision on the specific 
design and control of operating conditions. It was remarked that 
the connection between numerical modeling and experiments is 
too rare and also that material parameters in most cases are valid 
only for standard materials and not for the actual SOFC component 
microstructures. 

A review study on mathematical modeling based on steady- 
state and dynamic behavior, and control of polyelectrolyte 
membrane and SOFCs with respect to zero- one-, two-, and three- 
dimensional models was performed by M. Bavarian et al. [17]. The 
behavior of the fuel cell and processes that occur inside the fuel 
cells and contribute the existence of multiple time-scales in the fuel 
cells were investigated and also important components of these 
models were remarked. It was discussed that under which situa- 
tions a fuel cell observes steady state multiplicity and also stability 
of the steady state was issued. In continue, a review of control 
configurations and strategies by mentioning the advantages and 
disadvantages of each was presented. Finally, the topics that require 
further research studies are discussed. 

In this study, mathematical modeling of solid oxide fuel cells 
with respect to the tubular and planar configurations is reviewed 
into five subsystems considering the factors such as polarization 
losses, mass/energy/momentum conservations, diffusion through 
porous media, electrochemical phenomena in the PEN region and 
shift/reforming reactions inside them which is novel in this regard. 
Also, using variety of fuels fed to the SOFCs is discussed and their 
effect on the system is compared briefly. A short review of solid 
oxide fuel cell configurations and different flow manifolding are 
also presented in this study. For the first time in this area, novel 
models based on statistical data-driven approach existing in the 
literatures are studied. At the end of the paper, it is suggested and 
pointed which topics need more study for improving the fuel cell 
models. 

The remaining part of this article is organized as follows: Sec- 
tion 2 presents general principle of solid oxide fuel cell; Section 3 
discusses mathematical modeling of SOFCs, Section 4 addresses lit- 
erature on the SOFC configuration; and Conclusion is presented in 
Section 5. 


2. General principle of solid oxide fuel cell 
A solid oxide fuel cell uses a hard ceramic electrolyte and oper- 


ates at temperatures up to 1000°C. A mixture of zirconium oxide 
and calcium oxide forms a crystal lattice, although other oxide com- 


binations have also been used as electrolytes. The crystal lattice 
allows oxide ions to pass through it to reach the anode surface, 
where the oxide ions combine with H+ ions and form water. The 
solid electrolyte is coated on both sides with specialized porous 
electrode materials. The anode consists of a metallic nickel- and 
Y203-stabilized ZrO% skeleton, which inhibits sintering of the metal 
particles and provides a thermal expansion coefficient compara- 
ble to those of the other cell materials, thus limiting the buildup 
of stresses resulting from a difference in the coefficient of ther- 
mal expansion. The anode structure has a porosity of 20-40% to 
facilitate mass transport of the reactant and product gases. The 
cathode material mostly used is strontium-doped lanthanum mag- 
netite (La;_,SryMnO3, for x 0.10-0.15); its structure, such as that of 
the anode, is porous to permit rapid mass transport of the reac- 
tant and product gases. The cathode material has low levels of 
chemical reactivity with the electrolyte, which extends the life- 
time of the material. However, it is a poor ionic conductor, so 
the electrochemically active reactions are limited to the triple- 
phase boundary (TPB), where the electrolyte, gas, and electrode 
meet. In other words, the TPB of a fuel cell is the area of contact 
between the three phases necessary for electrochemical (hydro- 
gen oxidation and oxygen reduction) reactions at the electrode: 
an ion-conducting phase, an electron conducting phase, and a gas 
phase [18]. The bigger the TPB area of a cell, the better the qual- 
ity of the cell; a bigger TPB area allows the reactions to occur 
in more sites, thus maximizing current flow. Lanthanum stron- 
tium magnetite also works well as a cathode at high temperatures, 
but its performance quickly decreases as operating temperature 
decreases to less than 80°C. The electrochemical reaction occurs 
when both fuel and air are allowed to flow into the SOFC fuel cell 
separately. They diffuse through the porous electrode structure to 
the interlayer and are adsorbed. At the cathode interlayer, oxygen 
is reduced by the incoming electrons to produce oxygen anions 
that are conducted through the electrolyte to the anode inter- 
layer where they electrochemically combine with the adsorbed 
hydrogen to form water and release electrons to the external cir- 
cuit. 

The main advantages of the SOFC is that they can be operated at 
high efficiency of 50-60% and a separate reformer is not required 
to extract hydrogen from the fuel due to its internal reforming 
capability. Waste heat can be recycled to make additional elec- 
tricity by cogeneration operation [19,20]. The slow start up, high 
cost and intolerant to sulfur content of the fuel cell are some of its 
drawbacks. It is not suitable for larger fluctuations in load demand. 
Therefore, the SOFC is mainly used for medium and large power 
applications. In 1997 the Ceramic Fuel Cells Limited Company was 
demonstrated a5 kW laboratory prototype fuel cell system. Yakabe 
et al. [21] developed a 3 kW SOFC at Tokyo gas Co. Ltd. And they 
also analyzed the key factors to improve the performance of SOFC 
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Table 2 
The advantages and disadvantages of PSOFC and TSOFC. 


Type of SOFC Advantages Disadvantages 

PSOFC -Has the potential to offer higher power density than the tubular -The problems associated with the development and 
design. commercialization of planar SOFCs is requirement of high 

temperature gas seals. 

- Due to its compactness, it can be stacked in resemblance to polymer - Internal stresses in cell components due to non-uniform 
electrolyte membrane (PEM) fuel cells to satisfy the power temperature distributions and high manufacturing cost. 
requirement of an application. 
- Planar SOFC is simple to fabricate and can be manufactured into 
various configurations. 

TSOFC - Show great promise as potential candidates to replace conventional -Lower power density than the PSOFC design. 


heat engines. 
-Have higher energy efficiency. 


- Have high ohmic loss 


- They are designed for mid- to large-scale applications up to 2 MW. 


in the micro-grid system and are building a 250kW commercial 
stack model. 


2.1. SOFC configurations 


There are two principal SOFC configurations, tubular and pla- 
nar. The Schematic of tubular and planar SOFC is illustrated in 
Figs. 1 and 2 respectively. Tubular configuration consists of two 
tubes, an outer tube and an inner tube. The outer tube is the cell 
tube. The outer surface of the cell tube is anode side of the cell, 
and its inner surface is the cathode side. Between the anode and 
cathode sides lies the solid oxide electrolyte. The inner tube is an 
air injection and guidance tube, composed of alumina, from which 
preheated air is injected into the bottom of cell tube and flows over 
the cathode surface of the cell tube through the gap between the 
injection and cell tubes. The end of the cell tube is closed. Fuel gas 
flows over the anode surface through the gap among the cell tubes. 
Oxygen ions pass through the cathode and electrolyte and react 
with fuel, creating an electric current. A major advantage of tubu- 
lar Solid Oxide Fuel Cells (TSOFCs) over other types of fuel cells is 
that a variety of hydrocarbon-based gases or their synthesis deriva- 
tives, such as natural gas, biomass and coal can be potentially used 
as fuel sources. CO present in the synthesis gas can be oxidized in 
TSOFCs to generate CO% and electrical energy. The disadvantage of 
tubular SOFC compared to planar SOFC is that the ohmic losses are 
much higher due to the longer ways of electrical current. 

Generally, planar configuration is constructed of a positive- 
electrolyte-negative electrode (PEN), interconnect top and bottom 
of the cell stack, air and fuel channels. The fuel gas channel is 
between the anode and the separator plate and on the other side 
the air channel is located between the cathode and separator plate. 
Recognizing the advantages of the low-cost, simplicity in manu- 
facturing and high volume manufacturing with high volumetric 
power densities, most fuel cells manufacturers are concentrating 
on the planar SOFC concept. However, a challenge with the planar 
geometries is in obtaining mechanically stable structure, as thin 
layer ceramics are inherently susceptible to failure when subjected 
to moderate stresses. The planar concept provides the largest area 
for the current collector but it shows disadvantages in start-up time 
and mechanical stability. The advantages and disadvantages of two 
types of SOFC are depicted in Table 2. 

Based on the thicknesses of the electrodes and electrolyte 
layers, SOFCs are divided also into three types namely anode- 
supported, electrolyte-supported and cathode-supported. It was 
explored that ohmic and polarization losses in electrode-supported 
are smaller than in electrolyte-supported cells [22].However, the 
anodic polarization is smaller than the cathodic one [23]. More- 
over, anode-supported solid-oxide fuel-cells can operate at lower 
temperatures, which results in diminished manufacturing costs 
[24]. For these reason anode-supported cells are often considered 
to be the preferred configuration for SOFCs. But the thick anode 


Interconne 


Air 
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Anode (Ni-YSZ) 
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Fig. 1. Schematic of tubular SOFC. 


supported structure leads to enhance the concentration polariza- 
tion due to diffusion limitation, and thus decreases the fuel cell 
performance [25]. 

Regarding the flow configuration, there are three main types 
flow inside the fuel cells namely co-flow, counter-flow and cross- 
flow. 


3. Mathematical modeling of solid oxide fuel cells 


Mathematical modeling is an essential tool in designing fuel cell 
systems which represents the important aspects of the existing 
system and presents knowledge of that system in usable form. In 
other word, such models can provide a picture of voltage, current 
density, temperature, velocity and concentration of components 
as functions of position and time for various cell configurations 
and operating conditions. The effect of different variables and their 
sensitivity on the fuel cell performance can be studied by these 
mathematical models. A model is also useful in predicting the 
effects of altering process variables and in using that information 
to optimize cell performance. Moreover, dynamic simulation can 
provide reference for the control system research and design in 
future. 

Solid oxide fuel cell modeling for both the tubular and planar 
configurations tacked into account the electrochemical model, dif- 
fusion through porous media, the mass, energy, and momentum 
conservation balances. 

To perform first-principles model of the SOFC for the both con- 
figurations, single tubular and planar fuel cells are considered and 
divided into several subsystems as shown in Figs. 3 and 4. 

For the tubular SOFC the subsystems are as follow: 


e Subsystem 1 (SS1): air inside the injection tube; 
e Subsystem 2 (SS2): injection tube; 


S.A. Hajimolana et al. / Renewable and Sustainable Energy Reviews 15 (2011) 1893-1917 1897 


Interconne 
Anode (Ni-YSZ) 
Electrolyte (YSZ) 


Cathode (LSM) 
-aee à 
Air 


Interconne 


Fig. 2. Schematic of planar SOFC. 


e Subsystem 3 (SS3): air inside the space between the cell and 
injection tubes; 

e Subsystem 4 (SS4): cell tube; 

e Subsystem 5 (SS5): fuel mixture over the anode side of the cell 
tube. 


For the planar SOFC the subsystems are as follow: 


e Subsystem 1 (SS1): air channel; 

e Subsystem 2 (SS2): air side interconnector; 
e Subsystem 3 (SS3): fuel side interconnector; 
e Subsystem 4 (SS4): cell tube; 

e Subsystem 5 (SS5): fuel channel. 


The SOFCs models are then derived by writing the appropriate 
electrochemical and transfer phenomena equations. 


3.1. Model equations for subsystem 1 (SS1) 


This subsystem for TSOFC includes the air inside the injec- 
tion tube and for PSOFC is the air flow inside the air channel, 
between air side interconnector and cathode. In this subsystem, 
mass/energy/momentum transfer models can be accounted. 


3.1.1. Mass transfer model 

Gas transport is important to enhance the performance of the 
fuel cell. The performance of a SOFC at high operating current 
densities is reduced mainly due to the mass transport losses. 
Micro-structural parameters also largely affect the gas transport. 
Mass transport models inside the electrodes are giving an accurate 
prediction for gas concentration at the anode-electrolyte-cathode 
interface. Therefore; an accurate modeling of this phenomenon is 
imperative for development of better fuel cell designs. 


(i) A mass balance on the air inside SS1 for TSOFC takes the form 


api ie nt 

“at + Veg) = 0 (1) 
(ii) A mass balance on the air and O% inside SS1 (cathode channel) 

for PSOFC takes as the follow: 


apa” ch,c .ch,c 

air , > 

Ae + V(ugr Pair) =- VNo, (2) 
3p”: Yo 

ae V(U: p¥q, ) = —VNo, (3) 


where i denotes the generic i,, specie, Yo, the mass fraction of 
the O2 species, p is the density, u is the gas velocity and VNo, 
is the diffusion of O, into the porous cathode for participat- 
ing in electrochemical reaction, Superscripts inj and ch,c denote 
injection tube and cathode channel, respectively. 


Although mass transfer via diffusion inside the fuel cell sig- 
nificantly affects the fuel cell performance, some researchers like 
Costamagna et al. [26], Burt et al. [27], Stiller et al., [28], Ji et al. [29] 
neglected diffusion of components through media electrodes. Burt 
et al. [27] were assumed that changes in the x-direction are small 
and can safety be ignored. 

The diffusion through the porous media was modeled by writing 
detailed mass transfer equations with many researchers like Chan 
et al. [30], Yakabe et al. [31], Bove et al. [32], Ho et al. [33], Qi et al. 
[34] and others [35-41]. In these works, it was observed that the 
fuel cell behavior was more accurately than the models neglected 
the diffusion of gases in the electrodes. 

The diffusion of component through porous electrodes is 
assumed by using one of the three approaches as the following: 


(a) Fick’s Law model (FLM) is the simplest diffusion model and used 
in dilute or binary systems by Stiller et al. [28], Chan et al. [30], 
Bove et al. [32], Hussain et al. [37] and many other researchers 
[30,32,37-40,42-53]. The model is given by: 


N; = [Di V(pY;)] (4) 


where D; is the diffusion coefficient. 

(b) The Stefan—Maxwell model (SMM) is more commonly used in 
multi-component systems for calculation of specie fluxes. The 
model is given by: 


Y;N; — YiN; 
V(Y;) = ba i (5) 


where N is the species flux and C is the concentration. 

Some other researchers like Bhattacharyya et al. [54], Virkar 
et al. [55] Kim et al. [56] Bove et al. [57], Sanchez et al. [58,59] 
and Jin et al. [60] considered SMM model in their survey. 

(c) The Dusty Gas model (DGM), is the extended model of Stefan- 
Maxwell equation, which is the most complete diffusion model 
considering molecule-molecule and molecule-wall interactions 
in the porous media, commonly used as well, because it takes 
into account Knudsen diffusion phenomena and neglect the 
assumption of equal-molar counter diffusion in Fick’s law, 
which becomes invalid if the molecular weights of the diffusing 
gas species differ widely. The model is given by: 


iy , s YIN YIN; : 
=p 5+ Qo CD; (6) 


In many surveys this model is used [38,49,61-71]. Tseronis 
et al. [72] and Cayan et al. [40] used a combination of diffusion 
method. 

Tseronis et al. [72] discussed multidimensional model trans- 
port phenomena in an anode channel and diffusion through porous 
anode in an isothermal planer SOFC based on a combination of the 
Stefan-Maxwell and Dusty gas method. The results were validated 
with experimental data (Yakabe et al. [31]). 

Cayan et al. [40] illustrated a two-dimensional model based on 
mass transfer inside an SOFC. In this survey, diffusion methods, 
Stefan-Maxwell and Fick’s law relations, both including Knudsen 
diffusion, in order to predict the species concentration gradient was 
compared. The results showed that at low current densities, both 
models are well in agreement while by increasing the current the 
differences becomes bigger. 

A complete survey for performance comparison of Fick’s model, 
dusty gas model and Stephan-Maxwell model have been done by 
R. Suwanwarangkul et al. [52]. It was showed that DGM is the 
most suitable model for the H2-H20 and CO-CO2 systems due to 
the Knudson diffusion impact. However, due to its complexity, this 
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model can only be used at the conditions of low reactant concen- 
tration, high operating current density and small pore size where 
high accuracy of model prediction is needed. Otherwise, the SMM 
is a useful model for H2-H20 system and FM is suitable for the 
CO-CO2 system. SMM does not tacked into account the Knudson 
diffusion term in order to calculate the pore size effect while both 
FM and DGM do. It was fund that the equimolar counter diffusion 
is assumed to determine the flux ratio in both FM and the SMM. By 
contrast, the flux ratio in the DGM is calculated from the ratio of 
the square-root of the gas molecular weight. It was illustrated that 
DGM is suitable for multi-component system of Hyp-H20-CO-CO;, 
because the flux ratio calculation is complicated and equimolar 
counter diffusion cannot be assumed. Moreover, the Knudson dif- 
fusion effect must be considered. 

Diffusion through the porous material is typically described by 
either ordinary or Knudsen diffusion and has been found to play 
an important role in catalytic reaction. Ordinary diffusion occurs 
when the pore diameter of the material is large in comparison 
to the mean free path of the gas molecules. Molecular transport 
through pores which are small in comparison to the mean free 
path of the gas is regarded as a Knudsen type of diffusion. For 
Knudsen diffusion, molecules collide more frequently with the pore 
walls than with other molecules. Upon collision, the atoms are 
instantly adsorbed on to the surface and are then desorbed in a 
diffusive manner. As a result of frequent collisions with the wall 
of the pore, the transport of the molecule is impeded. The Knud- 
sen diffusion coefficient can be predicted using kinetic theory by 
relating the diameter of the pore and the mean free path of the gas. 
For straight and round pores, the diffusion coefficient is given by 
[30]: 


T, 
Dap = 97.07 , | — 7 
Ak Ma ( ) 


where Te is the cell tube temperature, and ř is mean pore radius 
which can be calculated from the bulk density, S, area of the 
porous solid, and ¢ porosity [30]. 


2e 


ro 
Sas 


(8) 

To account for the tortuous path of the molecule rather 
than along the radial direction and the porosity of the solid 
electrode for the fact that diffusion occurs only in the pore 
space, an effective Knudsen diffusion coefficient may be used, i.e. 
[30]: 


E 
Darceff) = Dak (3) (9) 


where ¢ is the tortuasity. 

Using the Chapman-Enskog theory of prediction, the binary 
ordinary diffusion coefficient in the gas phase [73] can be 
calculated as follows which was used by Yakabe et al. [21], 
Chan et al. [30], Ni et al. [41] and many other researchers 
[21,30,39-41,45,46,54,74-76]. 


(10) 


D 0.0018583 : 1) we 

- o.o018583( 4+ 1) 

mn Ma Mpg PO. Q DAB 
where p is the total pressure, 04g = (04 + Og)/2 is the collision diam- 
eter and pag is the collision integral based on the Lennard-Jones 
potential which can be obtained from eg, is the energy of molecular 
interaction: 


wo an 


Neufield et al. [73] is given by: 


1.06036 0.19300 
Be tiene " exp(0.47635(kT/€,p)) 
1.03587 1.76474 


* exp(1.52996(kT/éag)) | exp(3.89411(KT/€an)) (12) 


Hajimolana and Soroush [22], Qi et al. [34], Bhattacharyya et al. 
[54] and other authors [22,34,53,54,59,64,77-79] used Fuller et al. 
correlation which is given by [80]: 


1.013 x 10-7T.75[(1/Mo, ) + (1/Mn, )]'/ 


i= TPB cat cat 2 (13) 
[PIPE + PHC] — Co )ICKT/Eqp) 
The effective diffusion coefficient can be written as [80]: 
: eee (14) 


Dace) Dasteff)  Dakteff) 


Chapman and Enskog is the most common method for theo- 
retical estimation of gases penetration in the porous electrodes 
while Fuller correlation is derived based on the experimental and 
empirical data [81,82]. 

Table 3 shows some researchers that focused on mass transfer 
in their model. 


3.1.2. Heat transfer model 

In all the models that has been considered the thermal model, 
showed that the material properties, chemical kinetics, and 
transport properties of materials used in a single cell and stack 
depends on the temperature. Therefore, it is essential to develop 
a heat transfer model that could account for various heat effects 
in both of the solid structure and the flow channels. However, 
many researchers assumed an isothermal system in their model 
[38,60-62,64-66,83-87]. On the other hand, several literatures 
studied on heat transfer modeling in details in order to show the 
magnitude of heat transfer components. Fig. 5 shows a comparison 
between a non-isothermal and isothermal models. This figure 
depicts the responses of voltage in both isothermal and non- 
isothermal conditions by changing a fuel cell load. As it can be seen 
in Fig. 4, at the non-isothermal condition the voltage response is 
slow in comparison with the isothermal condition, since the heat 
capacity of solid materials is large and causing a slow response. 
This figure shows that temperature strikingly affects the fuel cell 
system and thus it is noteworthy to account energy balance in 
the SOFCs modeling. Table 4 depicts some authors that focused on 
heat transfer in their model. 


(i) An energy balance for the TSOFC inside the injection tube (SS1) 
has the form: 
o’ inj 
air air _ inj, inj pyinj inj inj 
Or!” VA Pairt air air + Yeon + rad (15) 
This equation represents the transport of energy in the injec- 


tion tube by bulk fluid flows (V A V ut Hi ), convection heat 
inj 


transfer (Qeony) and radiation between air flow components 
(Qah) 

ad a . . . 

Temperature of the inlet air stream into the SS1 has the 
strongest effect on the cell performance [22,43]. 

Convection heat transfer in SS1 for TSOFC is between air 
flow temperature and inside the wall of injection tube which 
is given by: 


Qidi = hin, (Tit, — TH) (16) 


where h is the convection heat transfer coefficient, T3 is the 
air flow temperature inside the injection tube. Subscript it; is 
inside the of injection tube. 
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Table 3 
Selected studies concerning on mass transfer through the SOFCs. 
Author Cell type Objective Method Short coming Assumption 
Izzo Jr et al. [38] Micro-TSOFC - Studying the gas flow Gas diffusion was The effect of DGM on fuel 1D, Steady state, 
through an modeled using the cell performance was not isothermal 
anode-supported tubular dusty-gas equations. compared with FL and 
SOFC and the subsequent SMM. 
diffusion of gas through 
its porous anode. 
-predict the polarization Mercury intrusion 
of the fuel cell porosimetry (MIP) was 
used to experimentally 
determine 
micro-structural 
parameters. 
Bhattacharyya TSOFC -Study the effect of In the process of -The effect of Three 2D diffusion, 
et al. [76] Knudson diffusion validation, phenomena diffusional models was isothermal dynamic 


Hussain et al. Anode-supported 


[37,83] SOFC 
Lehnert et al. [65] SOFC 
Suwanwarangkul TSOFC 

et al. [52] 


-Study the system 
response and compare 
with experimental 
results. 


-Develop a mathematical 
model to describe the 
transport of 
multi-component species 
inside the porous anode 
including reaction zone 
layer. 


-To anticipate the 
distribution of 
multi-component species 
in the electrode and 
reaction zone layers at 
different current density 
by the composition of 
ethanol reformate fuel. 
-Study the effect of shift 
reaction and finite 
reaction zone layer on 
the concentration 
overpotential. 

To describe the transport 
of gases inside the SOFC 
anode due to diffusion 


To evaluate the various 
mass transport models 
by validating them with 
the measured 
concentration 
overpotential data 
obtained from Yakabe 
et al. [42]. 


that affect the transient 
response of the cell 
significantly were 
identified. 


A modified 
Stefan-Maxwell 
equations incorporating 
Knudsen diffusion were 
used to model the 
multi-component 
diffusion inside the 
porous anode. 


Mass transport 
phenomena in the porous 
cermet, such as diffusion 
and permeation, 
accompanied by 
chemical and 
electrochemical 
reactions internal 
reforming and 
electrochemistry were 
studied and the resulting 
model is implemented in 
a numerical simulation 
program. 

Three types of models 
including Ficks’ model, 
the dusty gas model and 
the Stephan-Maxwell 
model inside the porous 
SOFC anode were 
developed to predict the 
concentration 
overpotential. 


not considered. 


-Further modifications 
may be needed in order 
to have a validated 
steady state and dynamic 
model in a broader 
operating range. 

-It was not clear among 
three diffuisional models 
why SMM was used in 
this model. 


-The momentum of 
gaseous species in 
porous anode was 
assumed to be negligible. 


-The effect of viscose and 
gravity for ideal gases 
was ignored. 


-It was not validated 
with experimental data. 


- The total pressure 
gradient inside the 
electrode was assumed 
to be negligible. 


model of an 
anode-supported 
tubular counter-flow 


1D steady state, 
Isothermal 


1D, isothermal steady 
state 


1D steady state, 
Isothermal 


(ii) An energy balance for the PSOFC inside the cathode channel 


(SS1) has the form: 
Of Hes ch,c, ch,c pych,c 
a S = VA pair Uair Hoir 


+ 


ch,c ch,c 


ad + Saif 


ch,c 
ony + 


(17) 


This equation represents the transport of energy in the cath- 


ode channel by bulk fluid flows (V Ap 


inj, inj ,yinj 
air Mair Hair 


), convection heat 


transfer (Qcony), radiation between air flow components, cathode 
side interconnector wall and cathode wall (Q,aq) and diffusional 
heat transfer by oxygen (Qgifr). 
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Fuel <4—— Fuel 
Cell tube 

Air Injection tube — ee 

Air ——» SS1 

Air SS1 
S82 ———— | Ss3 
SS4 
SS5 Fuel 

Fig. 3. Division of the single TSOFC into five subsystems. 
Anode side Interconnecto! SS3 
At the anode surface CH +H,O —} CO+3H) (Reforming) 
Fuel ———* SS5 Anode Electrolyte —— Depleted fuel CO+H,O 4P CO; H; (Shift) 
ss4 Cell Vv 
tube 5 Load 
Air ——> SSI Cathode —> — 
ir 
Cathode side Interconnector — [tH .., 
Fig. 4. Division of the single DIR-PSOFC into five subsystems. 
Table 4 


Selected studies concerning on heat transfer through Solid Oxide fuel. 


Author Cell type Objective Method Short coming Assumption software 
Inui et al. [89] PSOFC Numerically optimize the Concerning the calculation of The radiation heat transfer 1D transient 
operating parameters of the gas flow in each channel, was not considered 
the cell and the the one dimensional steady 
temperature of the inlet state mass and energy 
gas conservation equations were 
used, and the flow was 
determined by integrating 
these equations from the cell 
inlet to its outlet. 
Daun et al. [90] PSOFC Show the importance of Two different radiation solvers -The flow properties such 2D steady state 
radiation thermal transfer. were employed to solve for the as heat capacity, viscosity 
radiation source term in the were not considered, 
electrolyte: the Schuster 
Schwarzschild two-flux 
method [138] anda 
collision-based Monte Carlo 
method [139]. 
-The enthalpy and entropy 
were assumed constant 
and are not a function of 
temperature 
-Diffusional heat transfer 
was not considered 
Tanaka et al. [91] PSOFC Examine the effect of the It was developed a three The time dependent 3D steady state 
heat transfer due to dimensional simulation code of calculation was performed 
radiation on the cell the planar SOFC stack only for the cell 
operation, temperature distribution. 
Damm et al. [92] TSOFC Understanding, predicting, Modeling surface-to-surface Shift/ reforming heating 1D steady state 
and quantifying the effects radiative exchange, Radiative and diffusional heat 
of radiation in SOFC transfer in the flow Channels transfer were neglected. 
materials and systems. and Heat loss from the edges 
Brus et al. [93] TSOFC-IIR To describe numerical Heat transfer between the fuel -The flow properties such 
modeling of the radiative reformer surface and all other as heat capacity, viscosity 
heat transfer process surfaces facing the reformer were not considered 
surfaces was modeled 
-The enthalpies were 
assumed to be constant 
Sianchez et al. [58] TSOFC To simulate a complicated The work was divided in two -The velocity and mass 1D steady state 


heat transfer model and 
compare the errors used 
for those with simple 
models with the complex 
model. 


parts. First, two models for 
describing convective heat 
transfer was proposed. For the 
radiation heat transfer, two 
models were presented again 


flow were assumed 
constant along the length 
therefore the conduction 
heat transfer coefficient 
was constant. 
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Convection heat transfer in SS1 for PSOFC is between air flow 
and both cathode channel and cathode side interconnector walls 
which is given by: 


die = hint, (TSE — Tair) He het,co(Tct,co = Tair) (18) 


where subscripts int,c; and ct,Co are inside the cathode side inter- 
connector and outside the cathode side of the cell tube, respectively, 
and superscript cat is the cathode side. 

An empirical Nusselt number correlation is used to calculate the 
convection heat transfer coefficient: 


NU = f (Re, Pr) (19) 
hD 
NU = a (20) 


D is the fuel cell diameter. 
By Hayes and Kolaczkowski [88] expression Nusselt number 
yields as the follow: 


—0.5386 
NU = 3.095 4 8.933( E) ex (- ZA (21) 
Gz Gz 
where the Greatz number Gz is given by 
Gz= D Repr (22) 
x 


xis the axial position, Re and Pr are the Reynolds number and Prantl 
number, respectively. 

For a fully developed laminar flow approximation at constant 
wall temperature Nusslet number is assumed to be 3.66. 

For Re > 2300, 0.5 <Pr < 2000 and L> D, Gienelinski’s equation is 
applicable: 


Pr(Re — 1000)(f/8) Dn o) 
NU = ï4 23 
1+ 12.7(Pr2/3 — 1),/f/8 ( ( L ) ee 
f =(0.79 In(Re) — 1.64)! (24) 


where D,, is hydraulic diameter, L the length of fuel cell and ffriction 
factor. However this equation has complexity for calculating the 
convection coefficient. 

Colburn’s equations is easier to use but have errors around 20% 
[58]: 


For laminar flow: 
NU = 0.332Pr'/3Re!/2 (25) 
For turbulent flow: 


NU = 0.023Pr1/3Re4/5 (26) 


In some studies constant values for h and K at transient condi- 
tion are assumed [83,89]. In the study of D. S’anchez et al. [58] the 
effect of constant and variable convective heat transfer on fuel cell 
was compared. The results showed that these two models substan- 
tially affect the temperature at both ends of the fuel cell tube while 
there is no impact on the shape of temperature curve. Some authors 
like Hajimolana and Soroush [22], Ji et al. [29] and Qi et al. [43] cal- 
culated thermal conductivity, convection heat transfer coefficient, 
viscosity, specific heat capacity, enthalpy and specific resistivity as 
a function of temperature. The values of conduction and convection 
heat transfer coefficient existing in the literatures are illustrated in 
Table 5. 

Energy transportation inside and outside of the mass media 
occurs by diffusion of components. Diffusional heat transfer is 
available in [22,39,41,43,46,48,90-92]. However, this type of heat 
transfer was not involved in some literatures because it was 
assumed that it is not magnitude compared to the conduction heat 
transfer in porous media [61,67]. 


Table 5 
Different values of conduction/convection heat transfer coefficient. 
Authors he (kJ/m? sK) K (kJ/msK) 
Aguiar, P. et al. [95] hea =he,c=0.1 2x 10-3 
Damm et al. [92] hea =he,c=0.6 Ka =0.00584 
K-=0.00186 


Ketec = 0.00216 
Kair channei = 0.000047 


Kfuel channel =0.0002 
Denver et al. [92] hca = hec =0.302 


Neg =2.987 
hec = 1.322 


Xue et al. [84] Ka = Ke = Ketec = 0.00253 


For PSOFC in SS1, oxygen diffuses into the cathode in order to 
participate in electrochemical reaction, thus the diffusional heat 
transfer leads to: 


Qin = VNo, Ho, (27) 


Air mainly consists of non-polar nitrogen and oxygen molecules 
and therefore, is non-interaction (transparent) with radiation heat 
transfer at the moderate temperatures and pressures within SOFCs 
[93]. Therefore, no radiation heat transfer occurs between the air 
properties. So in SS1 for both TSOFC and PSOFC no radiation heat 
transfer occurs. 


3.1.3. Momentum model 

Velocity affects the film heat transfer coefficient on the anode 
and cathode side. An increase in the gas flow velocity leads to an 
increase in the film heat-transfer coefficient, causing the removal 
of more energy from the cell tube, resulting in a decrease in the 
cell-tube temperature [22]. On the other hand, velocity is affected 
by reactions, mass transfer, density change, etc. The flow velocities 
can be modeled by the momentum conservation. 


(i) Amomentum balance inside the TSOFC inside the injection tube 
(SS1) based on Navier-Stokes equation leads to [94]: 
inj , inj 
Ogi air) 
ot 
where ( V puu) is the momentum rate flowing in and out of the 
SS1, VP the pressure gradient, t is the shear stress, and F is 
prescribed body force (i.e. gravity force). 
(ii) A momentum balance for the PSOFC inside the cathode channel 
(SS1) based on Navier-Stokes equation leads to [94]: 


inj inj, inj, _ inj inj inj 
+ (V PGi Mair tair) SVE ON Ua + Fai (28) 


air air 


( ch,c chej 
air ~ air ch,c, ch,c, ch,cy __ ch,c ch,c 
or (V bait Uair Mair )= =No, Mo, Uair ~ V Pir 
ch,c ch,c 
E V Tair F Fair (29) 


where No, Mo, u is the momentum out flow rate resulting from 
O2 participation in the electrochemical reaction. 


Ferguson et al. [50], Janardhanan et al. [67] and other 
researchers [50,67,78,95-97] neglected the momentum con- 
servation in their model, while many researchers like Haji- 
molana and Soroush [22], Bove et al. [32], Autissier et al. 
[35], Qi et al. [43], Suwanwarangkul et al. [51] and others 
[32,35,39,43,45,51,53,54,61,63,66,69,77,78,83,86,91,98-103] con- 
sidered the momentum conservation in their model. Bhattacharyya 
et al. [54] compared the models with and without momentum con- 
servation and found that the change in velocity, especially in the 
cathode side channel is significant. When momentum conserva- 
tion is involved in the model the concentration of oxygen always 
remains higher, hence the oxygen concentration is also higher at 
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the TPB. It was observed that momentum conservation is important 
for high power range. 

The pressure gradient VP can be approximately determined 
using Poiseuille’s Law [104]. 

By calculating the friction factor f the wall shear stress yields: 


oo (30) 
pu 

fis constant for laminar flow and for turbulent flow it depends on 

Reynolds number and surface roughness [81]. 

In the case of the channel volume that is represented by porous 
medium (for the case of effective parameters) Darcy equation can 
describe the balance between the force from the pressure gradient 
VP and the frictional resistance € from the solid material which is 
given by: 
HgE 
K 
where Up is the pore velocity vector, € is the porosity (the volume 
fraction of void space), {tg is the dynamic viscosity of the gas mix- 
ture and K is the permeability which depends on micro-structure 
of the porous electrodes. 

For random packing systems of spherical particles, Caman- 
Kozeny correlation is generally used to evaluate the flow 
permeability as: 


-yP = Teg, (31) 


g3 


K = —— 
kg(1 — £) A 


(32) 
where the Kozeny constant kę is about 5 for porous media made of 
spherical particles. 

Ao is the specific area based on the solid volume: 


6 ne + (1 - Ney Jor? 


Ao = 
det Ney + (1 — Ne a3 


(33) 


a is bi-model size distribution. 

The volume flow rate per unit area through the porous structure 
is obtained by multiplying Up. The pore Reynolds number is defined 
by: 


ep = Ps |Up| dp (34) 


Hg 


where pg is the gas mixture density. Eq. (34) is valid for Rep < 1 [105] 
when pressure loss is dominated by pore friction and the inertia 
terms usually present in a momentum equation are insignificant. 
The Darcy equation describes the flow in the porous structure well 
away from the walls. It cannot model a no-slip condition at a wall 
nor the resulting boundary layers. It is found that, although the 
boundary layer thicknesses are typically only about 1-3% of the 
ceramic support thickness, Eq. (31) can be modified to allow their 
calculation. This is done by adding the so-called ‘Brinkman’ term to 
the Darcy equation which then becomes: 

HgE 
K 
where the normal stresses in the final term as usual are 
ignored. Darcy equation was used by several researchers like 
Yakabe et al. [21], Bove et al. [32], Ho et al. [46] and others 
[21,32,38,39,46,64,69,79,101-103,106,107]. 

Pressure drop in both anode and cathode sides was neglected 
by some authors like Shi et al. [85] Aguiar et al. [108,109] and 
others [78,85,96,97,108,110-113]. However, the temperature vari- 
ation along the cell can be accurately predicted by considering the 
pressure gradient among the fuel cell [114]. At the cathode side, 
neglecting the pressure gradient can lead to make an error in cal- 
culating of concentration overpotential by about 10% under typical 
working conditions [66]. 


-VP = ČE U, — V.(ugV)Up (35) 


3.2. Model equations for subsystem 2 (SS2) 


This subsystem for TSOFC includes the solid of injection tube 
and for PSOFC is the solid of the cathode side interconnector. In 
this subsystem mass and energy balances can be accounted. 


3.2.1. Heat transfer model 
(i) An energy balance for TSOFC in the injection tube (SS2) leads 
to: 
Cp.. Ti ; i i 
mi ge = onv + ond + Qaa (36) 
where m is the mass and Cp is the heat capacity. 

Convection heat transfer in SS2 for TSOFC is between inside 
injection tube and air flow and also between outside the injec- 
tion tube and the air flow inside the SS3 and can be written as 
follows: 


Qn = hit, (Toe — Tig) + hi, (Tair — Tito) (37) 


Heat transfers inside the gases and through the elec- 
trodes, electrolyte and interconnectors walls with conduction 
heat transfer, however, the thermal conductivity of gases is 
mainly very low (0.025 W/m K) [66]. Conduction heat transfer 
described by Fourier’s equation which is given by: 


Oop = VKYT (38) 


K is the thermal conductivity of the solid structure. 
Subscript it is the outside of the injection tube. 
Thermal radiation within tubular SOFC in SS2 occurs between 

injection tube and inside the cell tube walls which is given by 

eq (38): 

i o 

it = =— (Tå, = Ti) (38’) 

rad 
Raq is the radiation heat transfer resistance which is given by 

[80]: 

1-€¢ R 1 ; 1 — Eit Ait, 


Rrad = j u 
Ect Fet—it Eit Acti 


(39) 
(ii) An energy balance for PSOFC in the cathode side interconnector 
(SS2) leads to: 


3Cpint,c Tint,c 
Mint,c at 


= Qion + Qona + Qrad (40) 
Convection heat transfer for PSOFC is between the solid of the 
air side interconnector and air flow in SS2 which is given by Eq. 
(41): 
mie = hint,c;(Tair — Ta) (41) 
Thermal radiation within PSOFC in SS2 occurs between air 
side interconnector and cathode wall as follows: 


int, Oo 4 
git 2rd 5 oa (42) 
rad 
Ra 1 — Ecat , 1 = 1 = Eint Aint (43) 
Ecat Feat—int Eint Acat 


Because of the high operating temperatures of SOFCs (typically 
800-1200 K), there has been some estimation that thermal radi- 
ation may be an important factor of heat transfer within them. 
Radiation heat transfer has been involved by Hajimolana and 
Soroush [22], Damm et al. [89], Aguiar et al. [108], Iora et al. 
[91], Jia et al. [115], Achenbach [116] and other several mod- 
els in [26-28,35,44,51,61,75,77,83,89,91,96,97,108,110,115-120]. 
Detailed discussions of radiation heat transfer within SOFCs are 
given by Hajimolana and Soroush [22], Qi et al. [43], Ho et al. [46], 
Sánchez et al. [58], Damm et al. [89], Brus et al. [121], Wang et al. 
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[98], Calise et al. [122], Haynes [123]. However, some researchers 
ignored radiation heat transfer due to considering a uniform tem- 
perature everywhere in the cell [55], or because of assuming a thin 
wall of channels [124]. In Iwai, H., et al. [125] model it was not clear 
why the radiation heat transfer was not tacked into account. Wang 
et al. [126] and Yakabe et al. [31] neglected the effect of radiation 
heat transfer in their model due to a very small amount compared 
to other types of heat transfer. 

In fact, the magnitude of thermal radiation depends on two 
parameters, temperature and the thickness of electrodes and elec- 
trolyte. At a greater temperature and a thicker of the layers, the 
radiative transfer becomes more significant as compared to other 
heat transfer components. 

Calise et al. [122] presented a model of a tubular SOFC in which 
the radiation heat transfer was modeled in details. It was observed 
that radiation heat transfer has a strong effect on temperature dis- 
tribution among the fuel cell (about 70%). The model was developed 
by removing some simplification used in the study of Companari 
[127] and Stiller groups [28]. In Companari [127] the radiation 
heat transfer between air injection tube and the SOFC tube was 
neglected while F. Calise group [122] considered that the radia- 
tive heat transfer is predominately through the radial distribution 
among the fuel cell tube and air injection and is strikingly impor- 
tant. 

Sanchez et al. [58] performed a steady-state model of a tubular 
SOFC including convection and radiation heat transfers. Two mod- 
els for convection heat transfer were studied. In the first model 
film coefficient was a function of local temperature and compo- 
sition. The local value for the heat transfer coefficient and the 
thermal entry length of the cell is clear with this model. In the 
second model, an average value of the heat transfer coefficient 
for the whole length of the fuel cell channel was assumed. Based 
on this study, it was observed that a simple model can be used 
to evaluate the fuel cell performance with satisfactory accuracy, 
although there is a need to calculate local temperatures. For radi- 
ation heat transfer, two models were performed again. The first 
model accounted the radial radiation exclusively and thus, radiative 
exchange between adjacent cells is neglected. The second model 
considered the radiation completely in all directions. This causes 
to increase significantly the complexity of model. It was found that 
a simple model based on a constant film convective heat transfer 
coefficient and a simple radiation model is enough for predicting 
the fuel cell performance without internal information requiring 
and using a complicate model is not necessary. A complete heat 
transfer model, however, is necessary for evaluating the internal 
fuel cell temperature accurately. 

Conduction heat transfer was tacked into account by 
Stiller et al. [28], Negata et al. [44], Ho et al. [45,46], Fer- 
guson et al. [50], Sanchez et al. [58] and other researchers 
[28,44-46,50,58,63,67,75,77,78,89,92,99,101,102,108,110,115,117, 
128,129] in order to model the distribution temperature of fuel 
cell more accurately. However; no conduction heat transfer was 
used by Hajimolana and Soroush [22], Costamagna et al. [26], 
Burt et al. [27], Qi et al. [43] and Achenbach [116]. In this studies 
conduction heat transfer was neglected due to assuming a thin 
thickness of electrodes, electrolyte and interconnectors walls or to 
make simple the heat transfer modeling for calculating. It was also 
found that at low temperature the conductivity through the YSZ 
electrolyte becomes limiting [76]. 


3.3. Model equations for subsystem 3 (SS3) 
This subsystem for TSOFC includes the air inside the space 


between the cell and injection tubes and for PSOFC includes the 
fuel side interconnector. 


3.3.1. Mass transfer model 
(i) Amass balance on the air inside SS3 for TSOFC takes the form: 


apcat 

t+ Vien) =-VNo, (44) 
A mass balance on oxygen leads to: 

apet Yo 

eS + V(ustt oY, ) = —VNo, (45) 


The effect of O, diffusion on cell voltage and current is not sig- 
nificant because the concentration of O2 in the air is 21% [25,34]. 


3.3.2. Heat transfer model 
(i) An energy balance for TSOFC inside the SS3 leads to: 


cat pycat 


cat HC 
=a air = V Aptus eat a a A Qin (46) 


air “air “air iff 
which accounts for energy flowing into and out of SS3 with 
air(VA pat usat Heat), convectional heat transfer and diffusional 
heat transfer via oxygen into the cathode (Qi) 
Convective heat transfers from the cell and injection tubes 
are as follows: 
onv = Mito Titg — Tair ) + het; (Ter; — Tair ) (47) 


air air 


(ii) An energy balance for PSOFC inside the SS3 leads to: 


3Cpint,a Tint,a 
Mint,a OE 


inta , qint,a , qint,a , pinta 
=~ conv ! ond ' ad! ‘ond (48) 
Convection heat transfer for PSOFC is between the solid of the 
fuel side interconnector and fuel flow in SS7 and leads to: 
Qion = hint, a; (TR, = Te) (49) 


int,i fuel 


where superscript ano denotes the anode side. 
Thermal radiation within PSOFC in SS3 occurs between fuel 
side interconnector and anode side walls and is as follows: 
inta _ O 4 anot 
ad ~ R (Tet, T Titi ) (50) 
rad 
where subscript ct,4o denotes the outside the anode side of the 
cell tube. 


3.3.3. Momentum model 
Momentum conservation on the air inside SS3 for TSOFC takes 
the form: 
O( Pair uset) 
ot 


+ (V Pair stuctt) = —VNo, Mo, usat — y peat 


air air 


ay V reat a Feat (51) 


air air * 


3.4, Model equations for subsystem 4 (SS4) 


This subsystem for both TSOFC and PSOFC is the cell tube includ- 
ing electrodes and electrolyte. Electrochemical reaction, reforming 
and shift reactions occur in this subsystem. 


3.4.1. Electrochemical reaction 

High SOFC performance depends on optimal electrochemical 
reactions and mass transport processes. Since experimental stud- 
ies on SOFC are expensive, time-consuming and labor-intensive, 
quantitative mechanistic models for the cell PEN structure are 
essential for SOFC technology development. Electrochemical reac- 
tion occurs inside the fuel cell which as a result produces voltage 
and current. At cathode side oxygen ions (with a negative charge) 
migrate through the crystal lattice. The oxygen is supplied, usu- 
ally from air, at the cathode. Hydrogen and/or carbon monoxide 
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(CO) in the fuel stream reacts with oxide ions (O2~) from the 
electrolyte to produce water or CO and deposit electrons into 
the anode. The electrons pass outside the fuel cell, through the 
load, and back to the cathode, where oxygen from the air receives 
the electrons and is converted to oxide ions, which are then 
injected into the electrolyte. Generating efficiencies can range up 
to ~60%. In one configuration, the SOFC consists of an array of 
tubes. Another variation includes a more conventional stack of 
disks. The high operating temperature of the SOFC offers the pos- 
sibility of internal reforming. As in the molten carbon fuel cells, 
CO does not act as a poison and can be used directly as a fuel. 


H; +0% > H,0+2e 
Anode Side Reactions 4 CO + 027 —> CO; + 2e 
CH4 + 40°- > 2H20 + CO» + 8e 


Cathode Side Reaction {02(g) + 4e > 20%- 


The 1000°C operating temperature of SOFCs requires a signif- 
icant start-up time. The cell performance is very sensitive to the 
operating temperature. A 10% decrease in temperature can lead to 
a 12% drop in cell performance, because of a higher internal resis- 
tance to the flow of oxygen ions [46]. The high temperature also 
demands that the system include significant thermal shielding to 
protect personnel and to retain heat. Such requirements are accept- 
able in utility applications but not in small portable, transportable, 
or most transportation applications. 

It is significant that a SOFC can use CO, as well as hydrogen 
as its direct fuel. When both species of Hz and CO are present in 
the system, CO mostly participates in the water-gas shift reac- 
tion of equation [94] rather than in the electrochemical process 
[130]. However, current produces by the sum of hydrogen and car- 
bon monoxide oxidation at the anode side by the electrochemica 
reaction. The rate of H> oxidation is 2-3 times more than that of 
CO oxidation depending on the operationg temperature [131]. Ho 
et al. [132] assumed that the rate of CO oxidation is 3 times less 
that of H, oxidation due to high operating temperature of the fuel 
cell. On the other hand, Hofmann et al. [130] assumed that only 
H participate at electrochemical reaction due to its fast kinetic 
reaction. 

It is noteworthy that when the shift reaction is at equilib- 
rium condition, CO and Hz oxidation produces the same Nernst 
potential in the fuel cell [131]. Many authors like Cayan et al. 
[40], Yakabe et al. [42], Qi et al. [43] and other researchers in 
[40,42,43,61,68,69,112,116,133-135] considered CO oxidation in 
their model to make their model more realistic. However, some 
researchers like Costamagna et al. [26], Stiller et al. [28], Ho et al. 
[46] and also [26,28,46,102,103,106,110,113,117,120,136-138] 
ignored occurrence of CO oxidation at the anode side due to domi- 
nating of Hz over CO in the charge transfer chemistry It is estimated 
that about 98% of current is produced by H, oxidation in common 
situations and therefore CO seems to play a minor role and can 
safely be neglected [139]. 

Aloui and Halouani [69] compared the fuel cell performance 
with and without carbon monoxide oxidation at the anode side. It 
was observed by L. Petruzzi et al. [140] that CO oxidation at the elec- 
trochemical reaction influences the cell efficiency, fuel utilization, 
thermodynamic efficiency and electric power. Actually, oxidation 
temperature heavily depends on CO and as a result, CO oxidation 
can affect the cell performance. 

Methane can also possibly be electrochemically oxidized in the 
case of absence or near-absence of steam feed [93,104]. How- 
ever, with high content of steam the fast reforming reaction 
will dominate over any oxidation of methane and can safely be 
neglected. 


Incidentally, a pure voltage source is usually referred to as 
an electromotive force (emf). The electromotive force, reversible 
open-circuit cell voltage (denoted by Uocv) is given by the Nernst 
equation: 


0.5 
RT XH .aXQ, (Pe/ Pref) 
n 


Uocv = U%(T) 4 JF Te 
20,a 


(52) 


Equation [52] is hold for hydrogen, where Ugcy is the standard 
potential cell, R is the universal gas constant, F is Farady constant, 
P is the pressure, U°(T) is the standard cell potential, and x is the 
mole fraction. If the shift reaction is assumed to be at equilibrium 
the Uocv of Co is equal to that of CO. 

In practice, cell emf depends on temperature and concentra- 
tion of reactants and products. If the concentration of reactants 
increases relative to the products, the cell reaction becomes more 
spontaneous and the emf increases. As the cell operates, the reac- 
tants are used up as more product is formed causing the emf to 
decrease. 

For a fuel including Hz and CO the ideal potential for reaction 
using the Nernst equation is as follows: 


RT Py P, 
— y? De pery *CO 
Uocv = U'(T) + nF fns In ( P, ) +nco in ( P; ) 


P P P 
—nH, in( 2e) -no in | ran) (53) 


Ne is the total number of electrons transferred and given by [49]: 


Ne = Np, x 2+Nco x 2 (54) 


However, the actual cell voltage (U) is less than its theoretical 
open circuit voltage because it strongly affected by several irre- 
versible losses including activation losses due to irreversibility of 
electrochemical reactions at the three-phase boundary, concen- 
tration losses due to mass transport resistance in the electrode, 
especially for thick anodes as in an anode-supported SOFC [50] and 
ohmic losses due to ionic and electronic charge transfer resistances. 
Actual voltage is given by: 


U = Uocv — Nohm — Nact — Neonc (55) 


The rates of consumption of hydrogen and oxygen by the elec- 
trochemical reactions to generate an electric current of I are given 


by: 


LN; 
Rp, = (=): 
1\. 
Ro, = (7) 


The consumption of the reactants is accompanied by the pro- 
duction of water at the following rate: 


Rp, = (F) (57) 


3.4.2. Electrochemical losses 

3.4.2.1. Activation Losses. The activation polarizations are the 
results of the kinetics involved with the electrochemical reactions. 
It becomes an important loss when the current is low because at 
low current density (j) the reactants must overcome an energy 
barrier named activation energy (Esct) to drive the electrochemi- 
cal reactions at the electrodes-electrolyte interface and this barrier 
leads to the polarization. The activation barrier is the result of many 
complex electrochemical reaction steps where, typically, the rate- 
limiting step is responsible for the polarization. One can account 


(56) 
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for the anode and cathode activation polarizations using the well 
known Butler-Volmer equation [30]: 


nF afis nF 
j=io ex (Sa ) exo ( ( a )} (58) 


The equation explains the net anodic and cathodic current 
density due to an electrochemical reaction. Where subscript i rep- 
resents the species H2, CO and O3, j is the current density, n is the 
number of electrons transferred in the electrochemical reactions, 
Nact is the activation polarization, œ is the transfer coefficient and 
its value is usually 0.5, and jo is the exchange current. When £ =0.5, 
Eq. (58) can be expressed as follows: 


Fn a; 
j= 2jo sinh (= w) (59) 


2RT 
The total current density is given by: 


Jtot = jo, (60) 
Jtot = jn, +jco 


At the equilibrium potential the electrode/electrolyte interface 
still has electron transfer processes going on in both directions. 
The ongoing current in both directions at the dynamic equilibrium 
electric potential differences is called exchange current density. 
The exchange current density is a factor to measure the electro- 
catalytic activity of TPB for a given electrochemical reaction. It is 
also a significant factor for calculating the activation polarization. 
Exchange current density depends on many parameters such as 
concentration of reactants and products at the electrodes, temper- 
ature, pressure, the microstructure and electrocatalytic activity of 
the electrode, and even conductivity of electrolyte [104]. Exchange 
current density with respect to reactant concentrations depen- 
dency can be calculated as follows [30]: 


r o K 
Jo =Jo dl XK (61) 


where j? is a constant, and yx and Xx are the reaction order and mole 
fraction (at the electrode-electrolyte interface) for the Kth species, 
respectively [30,62]. A high exchange current density causes the 
reaction to occur rapidly, so a good fuel cell performance can be 
expected. 

Solving the Butler-Volmer equation with a transfer coefficient 
of 0.5 in the high-current-density regime (where the cell typically 
operates) leads to: 


BPS of Jf 
mar (E) fa) 


RT ees | J 
~(*_) sinh 
vee (nF) ae A Zoe 


(62) 


where joq and jo, are the anode and cathode exchange currents. 

Eact,ano and Eact,cat are energy activation of anode and cathode, 
PPB is the pressure of components inside the triple phase boundary, 
and A is the surface area of the cell. 

High activation polarization and low activation polarization are 
two possible cases of Eq. (58). Under high activation polarization, 
the first term of Eq. (58) will be much greater than the second term, 
so the simplified equation is given as follow: 


n= ($5) m() © 


which is called the Tafel equation. 
At low activation polarization, the term (Q@7qct,jF/RT) in Eq. (58) 
will be much less than and the exponential terms can be expanded 


as a Taylor series. Therefore, by ignoring the terms with order 
higher than unity the following equation yields: 


mni = (sa) J (64) 


niFjo 


which is a relation that shows linear current potential. The issue 
of applying the range of activation polarization for using the Tafel 
equation or linear current relation is discussed by S.H. Chan et al. 
[30]. 

Activation polarization in high temperature for the case of SOFCs 
is usually small [141]. It also depends on CO and H3 concentra- 
tion and decreases when the amount of these components in fuel 
increases due to the amelioration of the anode exchange current 
density. 


3.4.2.2. Concentration losses. Concentration losses are losses asso- 
ciated with concentration variation of the critical species due to 
mass transport processes. There are usually two sources of loss 
that are due to mass transport loss: (i) diffusion between the bulk 
flows and cell surfaces, and (ii) transport of reactants and prod- 
ucts through electrodes. Therefore, the concentration polarization 
is highly dependent on the gases used, as well as the distance 
through which the gases must diffuse. Pore volume percentage, 
as well as diffusion length, can be varied to optimize these prop- 
erties. For similar geometries, cathode concentrations are much 
larger than anode concentrations, because of the lower diffusivi- 
ties of O2/N2 in the cathode than H2/H20 in the anode. The anodic 
and cathodic concentration losses can be calculated, respectively, 
as the following: 


Nconc,a = neF 
es : (65) 

RT X; 
Nconc,c = >> (Vana,i — Vi) In (! il ) 


where X; is the component symbol for the ith species which may be 
involve in the fuel or oxidizer. v, and v; are the stoichiometric coef- 
ficients for the fuel and oxidizer mixtures, respectively, and Up ¡and 
vn, į are the stoichiometric coefficients of the ith product species 
in the fuel and air channels, respectively. [X;]* and [X;]’ denote 
the molar species concentration in the channels and triple phase 
boundary, respectively. 

This type of polarization loss is considerable for high fuel uti- 
lization [142] also normally reduced at high current density. A 
dramatic change happens for an anode-supported tubular SOFC 
which requires high fuel utilization due to a small error in the 
concentration overpotential calculation. 


3.4.2.3. Ohmic losses. Ohmic losses happen due to existence of 
resistance for ions transferring through the electrolyte and elec- 
trons through the electrodes and current collectors and by contact 
resistance between the cell components. Ohmic losses can be cal- 
culated as follows: 


Nonm = Roti (66) 


where Riot is the total area specific resistance and can be calculated 
as follows [142]: 


FORORO) o 


where ô is the thickness, and a, c and e denote the anode, cathode 
and electrolyte, respectively, and ø is the conductivity of electrodes 
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and electrolyte which is defined as a function of temperature as 
follows [143]: 


Oe = 3.34 x 104 exp (- a) 
z 95 x 10° 7 ( z 
a= T P(T 
42 x 10° 1200 
Oc = ( T ) exp ( T ) (68) 


The above equations are valid for the YSZ electrolyte material 
[144]. 

Hajimolana and Soroush [22] and Qi et al. [34,43] assumed 
a constant value for specific resistance while Costamagna et al. 
[26], Stiller et al. [28] and other many research studies calcu- 
lated the anode/electrolyte/cathode resistance [25,26,28,69,110, 
133,145]. 

In order to evaluate ohmic losses, some authors [28,59] used the 
Nisancioglu’s expression, Eqs. (69)-(74) [146], with a correction for 
the ionic/electric resistivity dependence on temperature, Eq. (74) 
giving by: 


A JjAsticeReq.j 


AVonm,j = ~F 5; 69) 
slice 
where 
Req = Req,1 “le Reg,2, 70) 
Rog (Pal Bal” + (P/S)? cosh(Ze) + PaPe/Sade(2 + Ze sinh(Ze))) 
eq, 1— j , 
2(1/pede)/*(pa/8a + pc/Se)"” sinh(Ze) oi 
Ree V Picmôicm(Pc/ôc) , (72) 
2 tanh(Licm/2/(Pc/Sc)/PicmSicm) 
a E S (73) 
e=3 Dede P at Pc/9c), 
T Ci 
j= — ex =}. 74 
Pi T P( T ) (74) 


In these equations A is the active area, j; the current density, 
Reqj is equivalent ohmic resistance, L is the slice length, ohmic 
resistance, ô thickness and C resistivity. 

Several SOFC models have been developed to study polarization 
losses within the PEN structure. 

Hirano et al. [146] presented a model in order to investigate 
the effect of concentration/ohmic polarization on electrochemical 
model by assuming that the activation polarization is negligible. 

Burt et al. [27], Ji et al. [29], Chan et al. [30], Izzo Jr 
et al. [38] and other many papers like [25,27,29,30,38,41,66, 
74,77,78,83,97,108,110,115,117,120,134,147-151] performed 
the three activatin/concentration/ohmic polarizations in their 
model. 

Tanaka et al. [152], Aguiar et al. [108] and Achenbach [116] 
considered only the activation polarization in their model. 

In the Aguiar et al. [108] survey, it was assumed that except 
for the shift reaction the principle electrochemical reactions are 
considered kinetically controlled, and therefore only the activation 
overpotentials were taken into account. It was assumed that the 
other overpotentials are small and can be ignored. 

Chaisantikulwat et al. [39], Ho et al. [45], Bhattacharyya et al. 
[54], Janardhanan et al. [67], Danilov et al. [99] and Wang et al. 


[98] accounted both activation and ohmic losses in their model and 
neglected the concentration polarization. 

Activatioin and ohmic polarization were accounted by Ota et al. 
[34]. In the study it was also assumed that concentarion gradient 
in the gas phase exists only in the porous electrodes. 

Padulleis et al. [153] neglected the activation/concentration 
polarization and only considered to ohmic loss in their model while 
Yakabe et al. [42] considered only concentration polarization in 
their model. 

In recent years, several SOFC models have been developed by 
Chan et al. [30], Sanchez et al. [59], Nam et al. [64], Klein et al. [90] 
and other studies in [30,55,59,62-64,85,86,90,136,148,150,154] to 
investigate reactions and transport phenomena within the PEN 
structure. These models differ widely in terms of their complexity 
and comprehensiveness. 

Chan et al. [30] performed a complete polarization model. In 
this paper the effect of electrodes and electrolyte thicknesses was 
studied. It was observed that the sensitivity of cell voltage due to 
change of electrolyte thickness is the highest, then sensitivity of cell 
voltage due to the change of cathode thickness and finally the sensi- 
tivity of cell voltage due to the change of anode thickness. Reducing 
all three component thicknesses causes, in general, reduced sensi- 
tivity strength with an improved operating current density range. 
However; considering both current density and PEN thicknesses 
is significant and needs an optimization between both factors in 
order to obtain a high performance while in this paper both fac- 
tors was issued. A cathode-supported cell is no way better than an 
anode-supported cell for high performance fuel cell design even 
under elevated pressure conditions at the cathode compartment 
to compensate for the serious cathode concentration overpotential 
loss. It was illustrated that the cell potential losses is mostly due to 
anode side of the cell in an anode-supported SOFC. 

Nam et al. [64] illustrated a comprehensive micro-scale model 
for explaining the transport phenomena and reaction in SOFCs. 
They showed that the potential loss in cathode was considerable at 
all current densities while the potential loss in anode was negligible 
at low current densities and important at higher current densities. 
Based on their research, it was observed that the size of particle 
has the most significant impact on the PEN performance optimiza- 
tion. By increasing the length of TPB in the electrodes and smaller 
particle size, the activation overpotential will be reduced, while at 
the same time by decreasing Knudsen diffusivity and flow perme- 
ability the mass transport resistance increases. On the other hand, 
the volume fraction of electronic phase of around 0.45 was found 
optimal for both anode and cathode for the same ionic and elec- 
tronic particle size. Higher porosity in anode and lower porosity 
in cathode gives a better PEN performance. The porosity, particle 
size, and thickness of electrodes were found to have conflicting 
effects on the PEN performance. For example, by decreasing the 
particle size of electrodes the efficiency of electrochemical reac- 
tion increases while the efficiency of mass transport decreases. 
An optimal set of micro-structural parameters, which resulted in 
best PEN performance, will not automatically result in best SOFC 
performance. Micro-structural optimization of SOFC performance 
requires consideration on the interconnect rib geometry and the 
flow field effects, in addition to the detailed micro-scale calcula- 
tion in PEN. In the paper also the effect of electrodes thickness was 
investigated in details. It was also depicted that with a larger of 
electrode thickness the performance increases, however; at higher 
current density the effect of electrode thickness will be negligi- 
ble. 

Jiang et al. [148] modeled a one-dimensional dynamic model 
of a TSOFC stack. It was found that higher operating temperature 
decreases both the Nernst potential and the irreversible losses, 
resulting in an initial increase then a decrease in cell efficiency. 
In this paper the resistivity and ohmic losses in different layers of 
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tubular SOFC were calculated as follows: 


p = 0.00298 e—1392/T fuel electrode 
p = 0.00294 e!9350/T electrolyte 


: 75 
p = 0.008114 e®°9/T air electrode (75) 
p = 0.12568 e4690/T interconnect 

L 


Jiang et al. also studied the effect of fuel cell size like fuel cell 
diameter and length. It was represented that by an increase of cell 
diameter the power increase due to a larger activation area at the 
same time and also the ohmic loss increase due to longer current 
path length. 

Ni et al. [66] developed an electrochemical model for studying 
the ammonia (NH3)-fed solid oxide fuel cells with proton- 
conducting electrolyte (SOFC-H) and oxygen ion-conducting 
electrolyte (SOFC-O). It was observed that the actual performance 
of the NH3-fed SOFC-H is significantly lower than SOFC-O, mainly 
due to higher ohmic overpotential of the SOFC-H electrolyte. It was 
also presented that compared with the NH3-fed SOFC-H, the SOFC- 
O has higher anode concentration overpotential and lower cathode 
concentration overpotential. The effects of temperature and elec- 
trode porosity on concentration overpotentials were also studied 
in order to identify possible methods for improvement of SOFC 
performance. This study observed that the use of different elec- 
trolytes not only causes different ion conduction characteristics at 
the electrolyte, but also significantly influences the concentration 
overpotentials at the electrodes. 

Ni et al. [41] investigated the effect of micro-structures of elec- 
trodes on SOFC performance. A mathematical model for a planar 
SOFC was developed for studying the effects of varies parameters on 
the fuel cell performance. Nam et al. [64] work, they revealed that 
both porosity grading and pore size grading can effectively improve 
the SOFC performance. Based on this survey, at low current density 
SOFC potential decreases with increasing porosity, while at high 
current density, optimal porosity was found and the optimal poros- 
ity increases with increasing current density. Moreover; the effect 
of pore size and porosity of electrodes on the overpotentials was 
studied while this effect was not considered by Nam et al. [64] work. 
Based on their results, variation of porosity has an inverse influence 
on activation and concentration overpotentials. Similar to poros- 
ity, the electrode pore size has opposite effects on activation and 
concentration overpotentials. The mechanisms of coupled trans- 
port and chemical reactions in the porous electrodes and evaluated 
possible ways to improve SOFC performance were also studied. The 
model developed in the paper can be effectively applied to the SOFC 
design optimization to maximize the energy efficiency of CH, fed 
SOFC for clean power generation. In addition; the impact of CH4 
fed SOFC on overpotentials was studied. Regarding this study, for 
SOFC fed by CHg, all the overpotentials decrease with increasing 
temperature while in Hj fed SOFC the concentration overpoten- 
tial slightly increases with increasing temperature. The effect of 
methane reforming and water-gas shift reactions is significant at 
high temperature, resulting in high rate of H2 production and high 
molar ratio of Hz and H20. 

Arpino et al. [74] investigated the effect of current density and 
fuel utilization on the overpotential losses. Based on their research, 
both concentration and activation polarizations at the anodic com- 
partment increase significantly for a high value of current density. 
It was also observed that there are no important gradients in the 
anodic concentration polarization at lower current density as the 
inlet hydrogen velocity varies. 

Table 6 illustrates some papers that concerned on PEN in their 
model. 


3.4.3. Mass transfer model 

At the electrodes, hydrogen and oxygen in the case of pure 
hydrogen as a fuel, diffuse from anode and cathode sides, respec- 
tively, into the triple phase boundary for occurring electrochemical 
reaction, and then after the products diffuse into the anode side. The 
diffusion of components from fuel channel into the TPB depends on 
the type of the fuel fed to the SOFCs. Mass transport models inside 
the electrodes are giving an accurate prediction for gas concen- 
tration at the anode-electrolyte-cathode interface. Therefore; an 
accurate modeling of this phenomenon is imperative for develop- 
ment of better fuel cell designs. 

A mass balance on oxygen inside the cathode-side diffusion 
layer yields: 


TPB 
Acat ə Po, Ro, 
R I | te) A un) 


Acat is the thickness of the cathode layer, R is the universal gas 
constant, P is the pressure of component at TPB, Ac is the surface 
area of cathode side and Tetypg is the cell tube temperature at TPB. 

By assuming that synthesis gas is as a fuel fed to SOFC, mass bal- 
ances on hydrogen and water vapor inside the anode-side diffusion 
layer yield are given as follows: 


PIPB 


Aano ð H2 Rp, 
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Aano is the thickness of the anode layer, Aa is the surface area of 
anode side. 


3.4.4. Heat transfer model 

Cell temperatures influence the electrochemical model and 
these temperatures affect the local driving voltage, polarizations, 
and heat generation within the cell. 


(i) An energy balance inside the SS4 for TSOFC leads to: 


Ipa Ta 
Met a T = Gav + Wg + Wy + We + Qer (80) 


Convection heat transfer inside the SS4 for TSOFC is between 
inside the wall of the cell tube and air flow inside the SS2 and 
outside the wall of cell tube and fuel flow which is given by: 


Qn = he; (TY — Te) + Neto (Thiel — Teto) (81) 


Thermal radiation within tubular SOFC in SS4 occurs between 
injection tube and inside the cell tube walls which is given by: 


o 
aa = pg Tito Tet) (82) 

For both TSOFC and PSOFC in SS4, O2 and Hz diffuse into 
the cathode and anode, respectively, in order to participate in 
electrochemical reaction, and H20 is the electrochemical pro- 
duction that diffuses from cell tube into the fuel channel. Thus 
the diffusional heat transfer leads to: 


TE 5 VN;H; (83) 


In the case of synthesis gas as a fuel, O2, H2, and CO diffuse 
into the cell tube, and H20 diffuses into the fuel channel. 

In the case of electrochemical heating (Qlec), heat gener- 
ates by the reversible and irreversible processes inside the 
electrodes and electrolyte within the triple phase boundary. 
Reversible heat generation is due to entropy transfer by the 
electrochemical reactions at the cathode and anode side which 


Xie and Xue. [87] 
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Table 6 
Selected studies concerning on PEN model of SOFCs. 

Author Cell type Objective Method Short coming Assumption 

F. Zhao et al. [158] PSOFC To anticipate the Porous media gas phase The models are limited to 1D 
concentration polarization transport processes based particular cells and cannot 
in detail. on Fick’s model was describe the effects of 

considered for predicting operating and design 
the concentration over parameters in detail. 
potentials. 

N. Autissier et al. [35] | PSOFC To predict current density, A commercial CFD tool was The model parameters are 1D steady state, 
flow, temperature and used, solving mass, not determined from laminar for both fuel 
concentration fields in momentum and energy experiments or from and air flow 
order to compare and equations; whereas literature, they were only 
optimize repeat element chemical kinetic equations estimated. 
geometry fora whole stack. were computed from 

external sub-routines, 
Radiative heat transfer was 
accounted between inner 
surfaces. 

Shi et al. [86] PSOFC To investigate the intricate The model was calibrated -Convection flux was 2D steady state, 
interdependency among using experimental neglected in the porous isothermal 


the ionic conduction, 
electronic conduction, 
multi-component species 
transport, electrochemical 
reaction processes and 
electrode microstructure 
for intermediate 
temperatures operation. 
between 750 and 850°C 


To de-convolute the 
complicated transient 
interactions among 
operating conditions, 
internal multi-physics 
processes, and SOFC 
current responses. 


anode-supported solid 
oxide fuel cell 


polarization curves and 
then validated by 
comparing each cell 
component polarizations 
(anodic, cathodic and 
electrolyte) 


The model was tacked into 
account the complex 
transport phenomena 
within PEN assembly, 
including charge 
(ion/electron) migration 
and species transport. 


electrode compared to 
diffusion. 


- Pressure gradients in the 
porous electrode were 
neglected. 

- All physical properties 
were evaluated at a fixed 
cell temperature. 

The pressure gradient 
within porous electrodes 
was neglected 


Isothermal 2-D 
transient 


is exothermic and slightly endothermic, respectively [155] 
which is given by: 


Qrev,a = Tetypp(2SH,0 — 2SH + 4Se-) Ap (84) 
Qrev,c = Tetzpg(2592- — 250, — 4Se-) ae (85) 


where Se is the molar entropies of the products and reactants 
at the operating temperature and pressure. 

Irreversible heat generation at the anode and electrolyte 
region is caused by anodic and cathodic overpotentials which 
are given by: 


Qirr,a = JNa (86) 


Qirr,c = jnc (87) 

Heat also generates by the ohmic heating in the electrodes 
and electrolyte due to existing of resistances to flow of ions 
and electrons which is given by: 


Qohm = Pej? (88) 


where pe is the ionic resistivity of the electrolyte. It is reported 
that ohmic heating is about 2.37-4.1% of the total heat pro- 
duced in an electrolyte-supported planar SOFC [156]. However, 
Andreassi et al. [129] assumed that the ohmic heating in the 
porous electrodes is negligible because of high electrical con- 
ductivity as compared to the ionic conductivity. 


All the heat sources generated by electrochemical reac- 
tion are determined by Bove et al. [32], Klein et al. 
[90], Ni et al. [103], Young et al. [104] and studies in 
[25,32,39,90,103,104,106,126,129,131,157,158]. 

(ii) An energy balance inside the SS4 for PSOFC leads to: 


OCpo, Tet 


t pt 


Me — 3t 


= onv 


| t | t | 
Oona F Qaa Qai + Qetc 


(89) 


Convection heat transfer inside the SS4 for PSOFC is between 
the wall of the cathode side and air flow inside the SS1 and the 
wall of the anode side and fuel flow which is given by: 


Oiv = Net,co(Tair = Tet,co) F hct,ao (Tarr = Tet,ao ) 


(90) 


Thermal radiation within planar SOFC inside the SS4 occurs 
between cathode side of cell tube and air side interconnector 
walls, and also between anode side of the cell tube and fuel side 
interconnector walls which is given by: 


ct _ 
ad Rrad 


on 


(romo 4 ) ; 
Rrad 


4 4 
int,i ct,do (Ter ) 


int,i ct,Co 


(91) 


In some works, radiation heat transfer was neglected in the 
planar SOFC because in these works it was shown that this type 
of heat transfer has only a negligible effect on the temperature 
field within the electrode and electrolyte [155,159,160]. 

Large variations in the cell performance occurs when only 
conductive and convection heat transfer among the fuel cell 
stack is considered, while by adding the radiation heat to the 
model, the uniformity of temperature promotes within the 
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Selected studies concerning on novel solid oxide fuel models. 


Author 


Objective 


Method 


Results and advantages of the model 


Entchey, E et al. [192] 


Wu et al. [197] 


Arriagada et al. [198] 


Yang et al. [199] 


To predict SOFC performance while 
supplying both heat and power to a 
residence. 


To describe the nonlinear 
temperature and voltage dynamic 
properties of the SOFC system 


To evaluating the SOFC 
performance 


To control fuel cell temperature 


Applied adaptive neuro-fuzzy 
inference system (ANFIS) 
techniques and artificial neural 
network (ANN) 


It was applied Takagi-Sugeno 
(T-S) fuzzy model. 


Used an Artificial neural network 
method 


A modified Takagi-Sugeno (T-S) 
fuzzy modeling study was 
presented. 


-The prediction of both ANN and ANFIS models 
was well in agreement with variety of 
experimental data. 


-Compared to the conventional models the 
ANN and ANFIS models are able to predict and 
optimize system performance faster and 
deliver better results in many instances. 

-the ANFIS model predicts the fuel cell current 
better than ANN model. 

-The ANN model predict the fuel cell voltage 
better than the ANFIS model. 


- The relationship between inputs and outputs 
of the SOFC system can be predicted while the 
complicated inner structure was ignored. 


-The SOFC model based on ANN was well in 
agreement with the physical model. 

-ANN model was much easier and faster to use 
compared to mathematical model. 


-The modified T-S fuzzy model simplified the 
complicated mathematical model and quickly 
made an accuracy relationship between inputs 


1909 


A radial basis function (RBF) 


and outputs of the system. 

-Based on the results the modified T-S fuzzy 
yielded a high accuracy over a relatively wide 
operating temperature range 


-The results were congruence in agreement 
with simulation results. 


genetic algorithm (GA) was 


applied. 


Wu et al. [200] To facilitate a valid control strategy 
design neural network based ona 
Hua et al. [201] -To facilitate a valid control 


strategy design. 


- Trying to avoid the complexities 
of model based on conversion laws 


Performed a black-box model of 
the SOFC based on least squares 
support vector machine 
(LS-SVM) 


- It was possible to design an online controller 
of a SOFC stack based on this GA-RBF neural 
network identification model. 


-The LS-SVM was superior to the conventional 
RBENN in predicting stack voltage with 
different fuel utilizations. 


stack, thus leads to more uniform the cell voltage. The varia- 
tion of cell-to-cell voltage is less than 0.2% in the model that is 
involved radiation heat transfer and 0.3% in the model without 
radiation [27]. In the case of indirect internal reforming, that 
radiative heat transfer accounts for up to 79% of the total heat 
transfer between the solid structure and the reformer [108]. 


3.4.5. Reforming and shift reactions 

Reforming reaction of hydrocarbons such as methane or natu- 
ral gas is a process for producing the hydrogen required by fuel 
cell. This reaction is predominately performed between 600 and 
900°C and is thus compatible with SOFCs. The heat produced in a 
fuel cell can provide the heat for the endothermic reforming reac- 
tion. Therefore, with elaborating an appropriate catalyst using in 
anode manufacturing, internal reforming reaction in fuel cell can 
eliminate the requirement for indirect internal fuel reformer and is 
expected to simplify the overall system design and finally the eco- 
nomical cost of these power systems will decrease substantially. By 
manipulating the internal reforming reaction the fuel cell tempera- 
ture can decrease to a suitable temperature. There are two alternate 
internal reforming within a SOFC: indirect (IIR-SOFC) and direct 
(DIR-SOFC) internal reforming. In the former process the reformer 
part is separate but adjacent to the anode cell. In the second process, 
methane is fed directly into the cell and the reaction occurs on the 
anode. However, this process has two major problems. One prob- 
lem is that internal reforming is highly endothermic and thus cell 


temperature drops initially and then increases again, hence gener- 
ate large temperature gradients across the cell, as well as a drop in 
cell performance. Deactivation of the anode by carbon formation at 
intermediate temperature and high pressure is the second problem 
of the DIR-SOFC which reduces the system efficiency. The problem, 
however, can be solved by adding H20 or CO3 into the syngas fuel. 
The complete steam reforming reaction takes place according to 
Eq. (92) 


CH4 + 2H, 0 => 4H C05 (92) 


Reforming reaction mainly occurs in the diffusional layer of 
anode or in the electrochemically active layer. But the rate of the 
reaction in the latter layer is assumed to be much lower than the for- 
mer layer because of competition with the electrochemical reaction 
which includes much more hydrogen. 

Eq. (92) is the combination of two independent reactions: 
methane steam reforming, Eq. (93), and water gas shifting, Eq. (96), 
which are highly endothermic and slightly exothermic, respec- 
tively: 


CH4 + H20 > 3H> + CO (93) 
CO + H20 > H; + CO2 (94) 


The shift reaction assumed to tacke place always at the anode 
layer. 
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Although Eqs. (93) and (94) occur, the following carbon forming 
reactions may also appear in the anode and the pre reformer: 


CH4 > C-+2H> (95) 


200 > C+ CO2 (96) 


These two reactions are endothermic and exothermic, respec- 
tively. If any of these two reactions take place at the anode, 
deposition of carbon particles over the anode surface will deac- 
tivate the catalyst and thus reduce the cell performance. At high 
temperature around 1000 K, methane is more likely to react with 
steam to form hydrogen according to Eq. (92) than to dissociate 
into atomic carbon and water [150]. This reaction increases when 
the steam concentration is high. For carbon monoxide, decomposi- 
tion occurs when the steam concentration is not high enough. Shift 
reaction is clearly dominant by the excess steam in the fuel stream. 

Many researchers like Hajimolana and Soroush [22], 
Ni et al. [41], Qi et al. [43], Nikooyeh et al. [48], Suwan- 
warangkul et al. [52] and also other researchers in 
[22,41,43,48,52,65,67,69,91,97,98,110,150,161,162] considered 
DIR-SOFC in their model. 

Yingwei et al. [97] observed that direct internal reforming pro- 
cess has significant effect on the distributions of the anode gas 
composition, cell temperature, solid temperature gradient, and cur- 
rent density along the cell length direction. 

Dokamaingam et al. [77,78], Brus et al. [121], Aguiar et al. 
[108,109] modeled IIR-SOFC in their system. 


3.4.5.1. Equilibrium model. Few authors considered that the 
reforming and shift reaction are equilibrium inside the cell 
[163,164]. The constant equilibrium for reforming and shift 
reactions yield with the [99,100] equations which depend on tem- 
perature exponentially as depicted in Eq. (100), where subscript i 
stands for both reforming and shift reactions. The constant values 
are reported in [163,164] which are given by: 


K Pi, Pco 
eq,ref = Dp 
Pcu,PH 0 (97) 
R _ Pm Pcos 
cast Pe oPco 
log Keg,i = a;T4 + biTZ. + CTS + diTct + e; (98) 


which both the reforming and shift reactions occur at the surface 
of the anode side. 

Suwanwarangkul et al. [61] calculated the shift reaction equi- 
librium constant as follows: 


Keg, shift = €XP(—0.2935Z? + 0.6351Z? + 4.1788Z + 0.3169) (99) 


where 
1000 
Z= TK) ~ 1 (100) 


However, equilibrium model is not a common assumption for 
refroming reaction in the papers. 


3.4.5.2. Kinetic model. Some researchers claim that equilibrium 
cannot be assumed for the reforming reaction inside the solid oxide 
fuel cell due to a slow reforming kinetics compared with both 
shifting and hydrogen electrochemical oxidation. Therefore, the 
following general expression is used to evaluate the reforming rate 
and depends on reactants concentration and temperature. 


E 
reng = APCH, * Pho XP (-a) (101) 


where E4 is the activation energy and À is the pre-exponential fac- 
tor. 


Hajimolana and Soroush [22], Qi et al. [43], Nehter et al. [75] 
and Aguiar et al. [108] assumed a kinetic reforming reaction in their 
model. 


3.5. Model equations for subsystem 5 (SS5) 


This subsystem for TSOFC is the fuel that occupies and flows in 
the space on the anode side of the cell tube and for PSOFC is the 
fuel channel. 


3.5.1. Mass transfer model 

The multi-component mass transport in SOFC anodes is mainly 
coupled with the bulk chemical reactions, diffusion and convection 
of species while at the cathode side there is no chemical reaction. 
A total mass balance on the fuel and component mass balances on 
the fuel components for both TSOFC and PSOFC leads to: 


ae? 
+ Vfl) = X VN: ue) 
dpan Yi 
Fe + (ugg Ppa Yi) = VN: + wi (103) 


w; is the rate of producing or consuming of components. 


3.5.2. Heat transfer model 
(i) An energy balance for the fuel gas inside SS5 for TSOFC results 
in: 
OP Fret H fuel 
ot 


In SS5, the convection heat transfer for TSOFC and PSOFC 
occurs between outside of the cell tube wall and the fuel flow 
inside the SS5 which is given by: 


= VA panguanorne + qano + Qane + Qao (104) 


soni = Neto (Tecto — Tiel) (105) 
Qa", in this subsystem for both TSOFC and PSOFC includes 
radiation heat transfer between walls and gas components. 
However, in most of the works, this type of radiation was 
neglected due to the lack of sufficiently accurate properties of 
gases. 

Fuel components including H20, CH4, CO, CO2 and possibly 
other hydrocarbons are very strange spectral, temperature and 
pressure dependence and so can be treated as transparent [93]. 

This heat flow per wall unit area is evaluated by using Eq. 
(106) where the emissivity of the mixture of gases £g is cal- 
culated as proposed by Hadvig in [165], Eq. (107) and its 
absorptivity ag as shown in Eq. (110) [166]. Subscripts w and 
g stands are for wall and gas in Eqs. (106)-(108) and i for the 
slice being analyzed. 


; Ew 4 4 
drad,i = 1-0- a1- al CeT gly i) (106) 
Eg = €CO2 — €H20 + €CO2EH20 (107) 
Tg 0.65 T, \ 0:45 
Og = eco, (= ) em0( z) £CO2£H20 (108) 
Tw Tw 


Few studies considered radiation heat transfer between walls 
and gas components in their model [59]. 
(ii) An energy balance for the fuel gas inside SS5 for PSOFC results 


in: 
appn H Ao 
ES = VA pfesuestiag +O ORR + Ome (109) 
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Voltage behavior at steady-state [83] 
Voltage behavior at transient state [30] 
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Fig. 5. The behavior of Fuel Cell Voltage at steady-state and transient condition. 


In SS5, the convection heat transfer for PSOFC occurs between 
both the anode side of the cell tube wall and fuel side intercon- 
nector and the fuel flow inside the SS5 which is given by: 


ano, = het ao(Tet,ao — Tet) + hint,q (T9? — Toe?) 


fuel fuel ( 11 0) 


The temperature distribution within the fuel cell is vital to 
SOFC operation, affecting the reaction kinetics and electro- 
chemical model and these temperatures influence the voltage, 
polarization, heat generation within the cell, thus the over- 
all cell performance. Based on the literatures, when only 
convection/conduction or convection/radiation heat transfer 
is considered large variations in cell temperature distribu- 
tion and then in the cell performance was observed. In 
Fig. 6, there is a comparison between three different mod- 
els regarding to distribution of cell temperature within the 
cell length. A complete heating components including radi- 
ation/convection/conduction/electrochemical and diffusional 
heating is involved in the T. Ota et al. [167] model, while R. 
Suwanwarangkul et al. [61] studied two different energy mod- 
els, first a model including only conduction, convection and 
radiation and in the second model the radiation heat transfer is 
ignored. It is observed that the temperature difference between 
the minimum temperature and maximum temperature within 
the cell length is 87 K in T. Ota et al. [167] model 05 K and 145 
in the first and second models of R. Suwanwarangkul et al. [61], 
respectively. Therefore, a strict temperature prediction distri- 
bution within SOFCs is essential for anticipating and optimizing 
the cell performances more accurately. 


3.5.3. Momentum model 
Momentum balance for the fuel in both TSOFC and PSOFC yields: 


Il Pfuel fuet) ano ,,ano,,ano ano ano 
at ae (V Pfuel fuel fuel) = 5 VN;MiU fiel = VPinel 


+ V Thiel + Fisel (111) 


Pressure gradient in the anode and cathode side is important to 
the electrochemical activities of the fuel cell. It was found that by 
ignoring the anode pressure gradient in the model can make and 
error in calculating of concentration overpotential by about 20% 
at a current density and temperature of 5000 A m~? and 1073K, 
respectively [66]. 


3.5.3.1. Different fuels fed SOFCs. One of the most brightened advan- 
tages of solid oxide fuel cells over the other types of fuel cells is that 
hydrogen and carbon dioxide are used as fuel in the cell. This means 


that SOFC can use many common hydrocarbon fuels such as natural 
gas, diesel, gasoline and alcohol without the need to reform the fuel 
into pure hydrogen. In other fuel cells, such as the polymer elec- 
trolyte fuel cell, which are fueled with pure hydrogen, the carbon 
dioxide is a poison. 

Hydrogen was considered as a fuel for SOFCs by Julian et al. [7], 
Wang et al. [126], Yuan et al. [168]. 

Hajimolana and Soroush [22], Autissier et al. [35], Aloui et al. 
[69], Aguiar et al. [108], Cui et al. [101], Janardhanan et al. [157], 
Fryda et al. [169] and Colpan et al. [170] considered hydrocarbon 
fuel fed to SOFCs. 

Many authors like Suwanwarangkul et al. [61], Aloui et al. [69], 
Xie et al. [86] and other researchers in [61,69,86,171-179] tacked 
into account biomass as a fuel for SOFCs because it is abundant and 
renewable. 

Hydrogen sulfide (H2S) is a byproduct of natural gas which can 
also use as a fuel in the solid oxide fuel cells. There are two possible 
overall reactions that use H2S directly in a SOFC: 


H2S + 502 <s H20+ js (112) 
3 
H25 + 502  H20 + S02 (113) 


However; it is not clear which of these reactions predominates 
in a SOFC running on H2S. 

Some researchers like Monder et al. [84], Liu et al. [163], Slavov 
et al. [164] and others in [84,163,164,180-185] used H2S as a fuel 
in SOFCs. 

Monder et al. [84] performed an isothermal steady state two- 
dimensional model of a planar SOFC fed by H2S. In this paper, the 
performance of fuel cell fed by H2S is not clear in comparison to 
other type of fuels. 

Ahmed et al. [161] were the first group that compared the effect 
of different fuels on SOFC performance. A cross-flow monolithic 
solid oxide fuel cell honeycomb (MSOFC) structure with alternat- 
ing layers of anode, electrolyte, cathode, and interconnect fed with 
different fuels, humidified hydrogen and coal gas, and air as oxi- 
dants was modeled. Results of the model with the simulated coal 
gas as the fuel and air as the oxidant showed that the MSOFC fed 
with coal gas reached a same cell voltage as a cell voltage of a hydro- 
gen fed MSOFC but the coal gas MSOFC was expected a lower net 
power and a lower average current density. 

Suwanwarangkul et al. [51], compared the performance of 
TSOFC using by biomass-derived synthesis gas as a fuel with the 
previous model (Suwanwarangkul et al. [61]) for a cell operating 
with humidified H2. The results depicted that the cell performance 
for operating with humidified Hz is more than that obtained with 
the synthesis gases because higher average temperatures within 
the cell were reached due to the higher energy content of humidi- 
fied H2. On the other hand, due to lower temperature distribution 
within the cell, the cell average temperature fed with synthesis gas 
is about 100°C lower than that of humidified hydrogen. 

Aloui et al. [69] presented an analytical modeling of polarization 
in a solid oxide fuel cell using biomass syngas product as fuel. The 
performance of fuel cell fed by biomass syngas was compared with 
the fuel cell fed by pure hydrogen. It was showed that a SOFC using 
the direct oxidation of syngas fuel (mixture of CO and H3) is more 
powerful than that using only H3 fuel. 

Niet al. [186] observed that the performance of CH, fed SOFC can 
promote at high operating temperature, as all the overpotentials 
decrease with increasing temperature. This behavior is different 
for Hz fed SOFC performance, because by increasing the fuel cell 
temperature the concentration overpotential also increases. 

P. Dokmaingam et al. [77], performed a mathematical two- 
dimensional steady state model of an indirect internal reforming 
tubular solid oxide fuel cell fed by four different primary fuels, i.e., 


1912 S.A. Hajimolana et al. / Renewable and Sustainable Energy Reviews 15 (2011) 1893-1917 


methane, biogas, methanol and ethanol. The objective of this paper 
was to investigate the effect of fuel type on the thermal coupling 
between internal endothermic reforming with exothermic electro- 
chemical reactions and system performance. It was observed that 
indirect internal reforming SOFC fed by methanol has the greatest 
performance in terms of electrical efficiency and temperature dis- 
tribution along the system. Based on this study, CO2 in biogas affects 
the system performance significantly. The electrical efficiency can 
be improved by removal of some carbon monoxide from biogas; 
however, a greater temperature gradient is expected. 

Brett et al. [187] were interested in studying on methanol fed 
SOFCs with an intermediate range of temperature between 500 and 
600°C. A lower operating temperature induces a more rapid start- 
up and shut-down, decreases the cell temperature fluctuation and 
then thermal stress in the ceramics. In addition, at a low operating 
temperature the dimensional thickness of the electrolyte can be 
reduced, giving almost the same results by promoting the ionic con- 
ductivity of the electrolyte at lower operating temperatures [159]. 
Moreover, fueling an SOFC directly by methanol ignores the inter- 
nal reforming process with large thermal gradients within the fuel 
cell. 

Arpornwichanop et al. [188] studied the performance of an 
anode-supported SOFC fed by ethanol under a direct internal 
reforming and isothermal condition. Electrochemical mode includ- 
ing all polarization losses was investigated in details. It was 
illustrated that cell temperature, the ratio of steam to ethanol and 
the inlet feed flow substantially affect the electrical efficiency of 
SOFC. In this study, the performance of fuel cell fed by ethanol was 
not compared with other fuels fed solid oxide fuel cell. 


4. Literature on the SOFC configuration 


Tubular and planer solid oxide fuel cells are discussed generally 
in Section 2. These two types of SOFCs are common configura- 
tions used in the literatures. Some researchers were interested in 
modeling both T-SOFC and P-SOFC. Ferguson et al. [50], studied a 
three-dimensional mathematical model on both tubular and plan- 
ner SOFCs with considering a local distribution of the electrical 
potential, temperature and concentration of the chemical species 
inside the both configurations. By comparing both tubular and pla- 
nar solid oxide fuel cells a less ohmic loss and higher efficiency for 
the former geometry was observed. 

Stiller et al. [28] like Ferguson’s group [50] were interested in 
both planar and tubular solid oxide fuel cells modeling. However; 
the study was focused more on effects of parameters on both fuel 
cell configurations with a 2D steady-state planar and tubular SOFCs 
fed by methane. The model was involved of ohmic resistance, con- 
vective, conductive and radiation heat transfers, chemical reactions 
and mass conservations. The required boundary conditions and 
input parameters of the models were discussed. The objective of 
the study was to observe the impact of inlet air flow temperature, 
pressure ratio, air flow rate and anode gas recycling on both con- 
figurations. It was illustrated that the cycle balance is different in 
both systems, as the tubular fuel cell required a lower air inlet 
temperature. Both fuel cell systems achieved above 65% electric 
efficiency. 

A complete comparison between tubular and planar solid oxide 
fuel cell configurations was presented by Zhang et al. [156]. A 
three-dimensional steady-state numerical model for both config- 
urations was described. This paper conducted detailed analyses of 
the electric polarizations for both configurations and presented the 
characteristics of the potential distribution in the electric com- 
ponents for the planar co-flow SOFC. The results were in a good 
agreement with the experimental data. Based on the results, for 
planar SOFC, polarizations results in the electrolyte ohmic loss and 


the activation polarization but for the tubular SOFC, the cathode 
ohmic loss has the greatest contribution to the sum of polariza- 
tions while the activation polarization is the next highest loss. It 
was also shown that the ohmic heat predominately depends on the 
current density in the planar SOFC while in tubular configuration 
the ohmic heat is produced from the cathode. The relative ohmic 
heat contribution to the total heat is higher for the tubular than 
for the planar configuration, and for the planar SOFC, this is more 
pronounced in the counter-flow design than in the co-flow design. 
Before Zhang et al. group [156] the effect of different flow manifold- 
ing, i.e., cross-co, or counter-flow were considered by Achenbach 
[116]. It was showed that the highest fuel cell efficiency is obtained 
for counter-flow, the most uniform current density distribution is 
reached for co-flow and the largest temperature gradients in the 
solid structure occur for cross-flow. Ferguson et al. [50] compared 
the counter flow and co-flow for planar geometry in their model 
and revealed that the counter flow for hydrogen as a fuel is the most 
suitable flow configuration. A comparison between counter-flow 
and co-flow was performed by Aguiar et al. [108]. It was presented 
that for the counter-flow configuration the temperature increases 
along the air flow direction with the maximum at the outlet and 
for the co-flow configuration the temperature profile reaches to 
its maximum value near the fuel flow outlet. In this system co- 
flow configuration is the more efficient one due to the high air 
temperature. 


5. Novel SOFC models 


As a kind of nonlinear multi-input—multi-output (MIMO) sys- 
tem, SOFC is hard to model by the traditional methodologies. 
Therefore; some researchers have attempted to establish novel 
fuel cell models by statistical data-driven approach. Data-driven 
approach is based on analyzing the data about a system, in partic- 
ular finding relationships between the input and output variables 
of the system without explicit knowledge of its physical behavior. 
The artificial neural network model is a novel model which is very 
useful for rapid calculation of SOFC e.g. dynamic, static and control 
simulations. The Hammerstein models is also a novel method for 
modeling SOFCs which include a nonlinear block followed by the 
linear system and are useful to model static and dynamic systems 
of the solid oxide fuel cell. 

Some researchers have attempted to establish novel SOFC mod- 
els [189,190,194-196]. 

Jurado [149] developed a multivariable Hammerstein model of 
a SOFC, which is applicable for small signal and dynamic stabil- 
ity studies. The Hammerstein models include a nonlinear block 
followed by the linear system and are useful to model static and 
dynamic systems of the solid oxide fuel cell. Dynamic behavior 
of temperature and output voltage losses was explained by con- 
sidering the electrochemical and thermal equations in this paper. 
Non-iterative algorithms for the identification of multivariable 
Hammerstein system were studied. The main debate in the deriva- 
tion of the results was the representation of the system using base 
functions. 

Bo Huo et al. [191] reported a dynamic modeling of SOFC using 
Hammerstein Model for dynamic simulation and control. The static 
nonlinear of the Hammerstein model was modeled by a radial basis 
function neural network (RBFNN), and the linear part was mod- 
eled by an autoregressive with exogenous input (ARX) model. In 
this study, a new gradient descent algorithm was used to antici- 
pate the hidden centers, the radial basis function widths and the 
connection weights of the RBFNN. For anticipating the parameters 
and the orders of the ARX model, the least squares (LS) algorithm 
and Akaike Information Criteria (AIC) were derived, respectively. By 
comparing the Hammerstein model and the RBFNN model exper- 
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Fig. 6. Comparison of cell temperature distribution between three different models. 


imentally, a better performance for the Hammerstein model was 
reported in this study. 

Wu et al. [192] presented a nonlinear offline model of SOFC 
used a radial basis function neural network based on genetic algo- 
rithm. The objective of the work was to optimize the parameters 
of radial basis function neural networks. The model was validated 
by the simulations. The simulation results illustrated that radial 
basis function neural network based on genetic algorithm method 
is superior to the conventional back progration neural network 
in predicting the stack voltage with different temperature basis 
parameter neural network in anticipating the fuel cell voltage 
in case of different temperature. Although many other operating 
parameters are also effective on fuel cell performance, only current 
density and temperature was considered in this model. Wu group 
[193] also studied the different flows effect on the performance 
of SOFC by developing an adaptive neural-fuzzy inference system 
(ANFIS). In order to identify linear and nonlinear parameters in the 
ANFIS, a hybrid learning algorithm combining back propagation 
and least squares estimate was adopted. The results were com- 
pared by simulations and the simulation results showed that the 
ANFIS model can efficiency predict the dynamic behavior of SOFC. 
However, the model was not included the temperature dynamic 
response and the SOFC fuel processor was also not considered. 

Milewski et al. [190] found that the Artificial Neural Network can 
be used to simulate the fuel cell behavior only by utilizing available 
experimental data. The cell voltage was predicted by the inputs of 
the model which were current density, temperature, fuel volume 
flow density, and oxidant volume flow density. The results reveal 
that the ANN can be efficiently used for the single SOFC with rela- 
tively high accuracy. Table 7 shows some researchers that focused 
on novel models. 


6. Conclusion 


Mathematical modeling is an essential tool in design of fuel cell 
systems, as it is important to understand the response of a cell stack 
under normal or transient conditions. Such models can provide a 
picture of cell stress, potential, current density, and temperature as 
functions of position and time for various cell configurations and 
operating conditions, and be used to examine the effects of changes 
on one or more variables and the relative system sensitivity to rel- 
evant design parameters. A model is also useful in predicting the 
effects of altering process variables and in using that information to 


optimize cell performance. Many research studies have been done 
in this regard. However, still many aspects of fuel cell modeling are 
remained for studying (Fig. 6). 

In this study, mathematical modeling of SOFCs for both tubu- 
lar and planar configurations was reviewed into five subsystems 
considering the factors like polarization losses, mass/energy/ 
momentum conservations, diffusion through porous media, and 
electrochemical phenomena in the PEN region and shift/reforming 
reactions inside them which is novel in this regard. Also, using vari- 
ety of fuels fed to the SOFCs was discussed and their effect on the 
system was compared briefly. A short review of solid oxide fuel cell 
configurations and different flow manifolding were also presented 
in this review study. For the first time in this area, novel models 
based on statistical data-driven approach existing in the literatures 
were studied. 

In this section we point to the subjects that need further 
research. 

Fuel cell temperature and voltage strikingly depend on electro- 
chemical reaction. The electrochemical oxidation of both hydrogen 
and carbon monoxide occurs on the triple phase boundary, whereas 
the proportion of these two reactions is still unknown. CO oxidation 
influences the overall performance of the cell and on the other hand, 
polarization resistances impress voltage directly. In many studies 
CO oxidation at electrochemical reaction is neglected or polariza- 
tion losses are not completely accounted. However, it is not still 
quite clear that by making some simplified assumptions for the 
electrochemical modeling how they can influence the accuracy of 
the model. What the difference between observed experimental 
values and the results of the mathematical modeling is. Calculating 
the errors even for small changes means that there is something 
we do not understand and need to account for it. Therefore, such 
phenomena should be deeply studied in both mathematical mod- 
eling and experimentation fields to perform a more accurate and 
validated model. 

Very limited studies focused on how the SOFC microstructure 
influence the chemical reactions and mass transfer in the porous 
SOFC electrode. Generally, length of TPB, porosity, particle size and 
thicknesses of electrodes and electrolytes has conflicting effect on 
the PEN performance. Thus, optimal design for PEN in order to 
improve the PEN efficiency is significant. 

Diffusion through porous media significantly affects the perfor- 
mance of solid oxide fuel cells. Fuel cell voltage and polarization 
losses depend on the concentration of components inside the triple 
phase boundary. In addition, diffusional heat transfer influences 
the fuel cell temperature. Therefore, it is necessary to consider 
diffusion in models in order to have an accurate prediction of 
fuel cell behaviors. Although several research has been done for 
presenting a beneficial diffusional model (FM, SMM, DGM) at differ- 
ent conditions, still more investigation requires to compare these 
three models together to minimize the error between mathemati- 
cal modeling and experimental results. 

Among the papers which were studied in this work, was shown 
that some of the models ignored radiative heat transfer between 
free walls and especial between gas components. Many works 
accounted radiation heat transfer in their model but the knowledge 
of thermal radiation property for both bulk and surface proper- 
ties at high temperatures relevant to SOFC is still incomplete and 
requires more study in this regard. 

Electrochemical heating is also still not fully available in open 
literature and is imperative to develop in the models. Ohmic heat- 
ing must be tacked into account in SOFC models due to the large 
ionic resistance of YSZ compared to heat of electrolyte materials 
used in other types of fuel cells. Moreover, conduction heat transfer 
in many papers was neglected while accurate determination of the 
temperature distribution including the radiation, conduction, diffu- 
sional heat transfer and electrochemical heating is required since 
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the material properties, chemical kinetics, and transport proper- 
ties of materials used in single cell and stack level depends on the 
temperature. 

Neglecting the momentum balance by several researchers was 
not fully justified. More research studies require showing the effect 
of momentum balance especially inside the porous electrodes on 
solid oxide fuel cell performance. 

The pressure gradient in the electrodes of SOFCs is simply 
ignored without any justification in modeling studies. Therefore 
the effect of pressure distribution on mass transfer and thus on 
overpotential should be more studied. 

A major advantage of solid oxide fuel cells over other types 
of fuel cells is that a variety of hydrocarbon-based gases or their 
synthesis derivatives, such as natural gas, biomass and coal can 
potentially be used as fuel sources. It was shown that a SOFC using 
the direct oxidation of syngas fuel (mixture of CO and H3) is more 
powerful than using pure H3 fuel, and also methanol in compar- 
ison to other fuels like methane, biogas, and ethanol as fuels, has 
the greatest performance in terms of electrical efficiency. However, 
still a complete research study needs to present a suitable fuel for 
promoting the SOFCs performance and their power efficiency. In 
addition, the effect of each fuel on the important factors of SOFCs 
such as cell temperature, voltage, reactions, and cell performance 
is not completely studied. 

In many papers the effect of reformers and different geome- 
tries and configurations is not clear. It is imperative to show the 
importance of the mentioned factors in the models. 

The intelligent novel models avoid the complexity of the math- 
ematical model and provide the input-output relationship much 
faster and easier. Based on these models, valid control strategy 
studies like predictive control and robust control can be developed. 
Therefore, more research studies requires for developing the novel 
models. 
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